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ABSTRACT 



We present new synthetic models of the TP-AGB evolution. They are computed for 7 values of initial metal content (Z from 0.0001 
to 0.03) and for initial masses between 0.5 and 5.0 M Q , thus extending the low- and intermediate-mass tracks of Girardi et al. (2000) 
until the beginning of the post-AGB phase. The calculations are performed by means of a synthetic code that incorporates many 
recent improvements, among which we mention: (1) the use of detailed and revised analytical relations to describe the evolution of 
quiescent luminosity, inter-pulse period, third dredge-up, hot bottom burning, pulse cycle luminosity variations, etc.; (2) the use of 
variable molecular opacities - i.e. opacities consistent with the changing photospheric chemical composition - in the integration of a 
complete envelope model, instead of the standard choice of scaled-solar opacities; (3) the use of formalisms for the mass-loss rates 
derived from pulsating dust-driven wind models of C- and O-rich AGB stars; and (4) the switching of pulsation modes between the 
first overtone and the fundamental one along the evolution, which has consequences in terms of the history of mass loss. It follows 
that, in addition to the time evolution on the HR diagram, the new models predict in a consistent fashion also variations in surface 
chemical compositions, pulsation modes and periods, and mass-loss rates. The onset and efficiency of the third dredge-up process are 
calibrated in order to reproduce basic observables like the carbon star luminosity functions in the Magellanic Clouds, and TP-AGB 
lifetimes (star counts) in Magellanic Cloud clusters. In this paper, we describe in detail the model ingredients, basic properties, and 
calibration. Particular emphasis is put in illustrating the effects of using variable molecular opacities. Forthcoming papers will present 
the theoretical isochrones and chemical yields derived from these tracks, and additional tests performed with the aid of a complete 
population synthesis code. 

Key words. Keywords should be given 



1. Introduction 

Owing to their large luminosities and cool photospheres, ther- 
mally pulsing asymptotic giant branch (TP-AGB) stars are 
among the most remarkable objects in near- and mid-infrared 
surveys of nearby galaxies. The TP-AGB is also one of the most 
intriguing phases of stellar evolution, marked by the develop- 
ment of high mass loss rates, the presence of recurrent lumi- 
nosity and temperature excursions, the rich nucleosynthesis, and 
the long-period variability. Several "third dredge-up" episodes 
- followed by hot-bottom burning (HBB) in the most massive 
AGB stars - cause the surface pollution with He and CNO ele- 
ments of primary origin, which become directly observable first 
via the changes in the spectral features of AGB stars themselves, 
and later in the emission line spectrum of the subsequent plane- 
tary nebulae phase. In the context of galaxy models, AGB stars 
are crucial contributors to both their integrated spectra and their 
chemical enrichment. 

Computations of the TP-AGB phase by means of complete 
evolutionary codes are very demanding in terms of computa- 
tional time, and in most cases they fail to predict basic ob- 
servational facts such as the conversion from M to C spectral 
types in AGB stars with initial masses and luminosities as low 

as Mi ~ 1.5 M and M\, Q \ 3.5 (see Herwig 2005 for a review, 

and Stancliffe et al. 2005 for a recent exception). Such difficul- 
ties derive from the complex structure of these stars, as well as 



the uncertainties in the modelling of convective dredge-up and 
mass-loss processes. In order to provide extended grids of TP- 
AGB models that reproduce basic observational constraints, the 
only viable alternative so far is the use of the so-called synthetic 
codes, in which the stellar evolution is described by means of 
simplified relations derived from complete stellar models, while 
convective dredge-up and mass-loss are tuned by means of a few 
adjustable parameters. 

Purpose of this paper is to describe new improved models of 
the TP-AGB phase of low- and intermediate mass stars, whose 
computation has been motivated by a series of recent improve- 
ments in the area: 

1) Starting from the work by Wagenhuber & Groenewegen 
(1998), recent theoretical works (Karakas et al. 2002, hereafter 
K02; Izzard et al. 2004, hereafter 104) have significantly detailed 
and extended the range of validity of analytical relations de- 
scribing the evolution of TP-AGB stars - e.g., the core-mass 
luminosity relation, inter-pulse period, pulse cycle variations, 
dredge-up efficiency, etc. The theoretical modelling of the mass- 
loss phenomenon has also significantly improved (Willson 2000; 
Winters et al. 2003; Hofner et al. 2003, and references therein). 
Moreover, Marigo (2002) has recently demonstrated that the 
low-temperature opacities used in present-day TP-AGB calcu- 
lations are in substantial error, due to the non-consideration of 
the dramatic changes in molecular concentrations that occur in 
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the outer envelope of TP-AGB stars as their chemical compo- 
sition is altered by third dredge-up episodes; the consequences 
of using consistent molecular opacities turn-out to be remark- 
able when stars enter into the C-rich phase of their evolution 
(see Marigo 2002; Marigo et al. 2003). All these theoretical im- 
provements should allow a better description of many aspects of 
the TP-AGB evolution, and also a better understanding of their 
dependence on the stellar mass and initial chemical composition. 

2) Present-day near-infrared cameras are enormously in- 
creasing our knowledge of the AGB populations of Local Group 
galaxies. DENIS and 2MASS have provided photometry in the 
red and near-infrared for the complete sample of non-obscured 
AGB stars in the Magellanic Clouds, revealing striking features 
such as the "red tail" of carbon stars (Cioni et al. 1999; Nikolaev 
& Weinberg 2000). Dedicated C star surveys (e.g. Battinelli & 
Demers 2005ab, and references therein), as well as future near- 
IR surveys using UKIDSS and VISTA, will soon provide com- 
plete TP-AGB samples for a large subset of the Local Group. 
On the other hand, infrared surveys (IRAS, ISO, Spitzer) are re- 
vealing the samples of dust-enshrouded, optically obscured, TP- 
AGB stars in the Magellanic Clouds (Loup et al. 1999; Cioni et 
al. 2003; van Loon et al. 2005; Blum et al. 2006). Interpreting 
these large photometric databases in terms of galaxy properties 
(star formation rate, age-metallicity relation, density profiles) re- 
quires the use of suitable TP-AGB models. An example of the 
potentialities of present-day data-sets are given by Cioni et al. 
(2006ab), who detect variations of the mean stellar metallicity 
and star formation rate across the LMC and SMC galaxies, us- 
ing only the near-IR properties of their AGB stars as compared 
to those predicted by theoretical models. 

3) Micro-lensing surveys in the Magellanic Clouds and in 
the Galactic Bulge (MACHO, OGLE) have provided a huge 
amount of high-quality optical data for long-period variables; 
these data are nicely complemented with near-IR photometry 
from DENIS and 2MASS. Particularly striking has been the dis- 
covery of several different sequences in the luminosity-period 
plan, four of which represent the various pulsation modes of 
AGB and upper-RGB stars (Wood et al. 1999; Soszynski et al. 
2006; Groenewegen 2004; Ita et al. 2004; Fraser et al. 2005). 
The data clearly indicate that TP-AGB stars start pulsating as 
small-amplitude semi-regular variables (first and second over- 
tone, and mixed modes) and then later in the evolution become 
high-amplitude Mira variables (fundamental mode). Since it is 
generally believed that the pulsation period correlates with the 
mass loss rate of such stars, this may have important evolution- 
ary implications that still have to be explored. 

Therefore, we aim at producing TP-AGB models - and 
their derivatives, like isochrones, luminosity functions, syn- 
thetic colour-magnitude diagrams, etc. - including the above- 
mentioned improvements in the analytical relations, molecu- 
lar opacities, and mass-loss prescription, and consistently pre- 
dicting variations in surface chemical compositions, pulsation 
modes and periods, and mass-loss rates. Like in any other set 
of synthetic models, uncertain parameters are to be calibrated 
in order to reproduce basic observables like the carbon star lu- 
minosity functions in the Magellanic Clouds, and TP-AGB life- 
times (star counts) in star clusters. In a second phase, we intend 
to check also relatively new observables like the period distri- 
butions and relative numbers of stars in the several sequences 
of LPVs, and the properties of dust-enshrouded AGB stars with 
high mass-loss. 

The present paper is organised as follows. In Sect. [2] we de- 
scribe model ingredients and computations. Sect. [3] illustrates 
the basic model predictions, with special emphasis put on the 



distinctive features found in the present calculations. Section [4] 
together with the Appendix detail the calibration of model pa- 
rameters, based on selected observational constraints. Finally, in 
Sect. [5] we briefly recall the major achievements of this work, 
and the next steps to be followed. 

2. Synthetic AGB evolution 

A very detailed code for the synthetic evolution of TP-AGB 
stars has been developed by Marigo et al. (1996, 1998) and 
Marigo (2001 and references therein). It couples the use of up- 
dated analytical relations (e.g., the core mass-luminosity rela- 
tion, the interpulse duration), with a parametric description of 
the third dredge-up episodes, and numerical integrations of a 
complete envelope model. Over the last years, several additional 
improvements have been introduced, the most relevant of which 
are briefly recalled: 

- Development of a consistent method for dealing with the 
over-luminosity effect caused by hydrogen burning at the 
bottom of the convective envelope and its related nucleosyn- 
thesis (HBB; see Marigo 1998). 

- Adoption of more physically-sound dredge-up parameters, 
linked to the minimum post-flash temperature at the bottom 
of the convective envelope, T dred . The two free parameters 
associated with the possible occurrence - described by T dred 
-, and the efficiency - described by A - of the third dredge-up 
were calibrated by fitting the carbon stars luminosity func- 
tions in both Magellanic Clouds (Marigo et al. 1999). 

- Calculation of molecular opacities properly coupled to the 
actual chemical composition of the envelope (A\, ar ), in place 
of the common and inappropriate choice of low-temperature 
opacity tables (/Cfi x ) strictly valid for solar-scaled chemical 
mixtures (Marigo 2002; Marigo et al. 2003). 

The present study stands on these previous works, adopting 
additional improvements to be described below. For a general 
description on how synthetic TP-AGB models work, the reader 
is referred to Groenewegen & Marigo (2003). 

2. 1 . Initial conditions 

The initial conditions at the first thermal pulse - namely total 
mass, luminosity, core mass, and surface chemical composition 
as affected by the first and second dredge-ups - are taken from 
Girardi et al. (2000) for the initial metallicities Z = 0.0004, 
0.001, 0.004, 0.008, 0.019, and 0.03. For Z = 0.0001, we use 
an additional set of tracks computed by Girardi (2001, unpub- 
lished; see http://pleiadi.oapd.inaf.it) using the same 
input physics as Girardi et al. (2000). Without entering in de- 
tails, suffice it to recall that these tracks are computed with OPAL 
(Iglesias & Rogers 1993, and references therein) and Alexander 
& Ferguson (1994) opacities, and include a moderate amount of 
convective core overshooting As a consequence of overshoot- 
ing, the upper mass limit for the development of the AGB phase 
is located close to 5 M Q (instead of the ~ 8 M found in classical 
models) for all values of Z, and is slightly uncertain because the 
initial phase of carbon burning has not been followed in detail 



1 Namely, in these models the overshooting parameter A c (see 
Alongi et al. 1993) is set to 0.5 for initial masses M/ > 1.5 M e , which 
implies an overshooting region extending for about 0.25 pressure scale 
heigths above the classical Schwarzschild border. For M/ < 1.5 M Q , A c 
is gradually reduced until it becomes null for Mj < 1.0 M Q 
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(see section 3.1 in Girardi et al. 2000). Therefore, Mi — 5 M Q is 
the largest initial mass to be considered in this work. 

Finally, we remark that throughout this paper, except when 
otherwise stated, the term initial mass is equivalent to the stel- 
lar mass at the first thermal pulse, as the pre-AGB evolution in 
Girardi et al. (2000) is followed at constant mass. The effect of 
mass loss on lower-mass models while evolving on the red gi- 
ant branch will be considered a posteriori in the construction of 
stellar isochrones (Marigo & Girardi, in prep.). 

2.2. Luminosity and core-mass growth 

Important improvements have been introduced in the present 
calculations to follow (i) the evolution of the stellar luminos- 
ity L(t) and, (ii) the growth rate of the core mass dM c /dt = 
Lh Qh/X(H), where Lu is the rate of energy generation due to 
radiative shell H-burning, Qu is the emitted energy per unit mass 
of H burnt into He, and X(H) is the H abundance in the stellar 
envelope. 

To this aim we need first to choose a core mass-luminosity 
(M c - L) relation, among the various available in the literature. 
This key-ingredient of any synthetic AGB model is usually de- 
rived from complete AGB evolutionary calculations by fitting 
the surface luminosity at the stage of the quiescent pre-flash 
maximum as a function of the core mass, and possibly including 
the effect of other less influent parameters (metallicity, envelope 
mass, initial conditions at the first thermal pulse). 

For lower-mass models not experiencing HBB, we adopt 
the core mass-luminosity (M c - L) relationship provided by 
Wagenhuber & Groenewegen (1988; their equation 5), that in- 
cludes important effects like dependence on metallicity and the 
peculiar behaviour owing to the first thermal pulses. 

The deviation from the M c - L relation due to the overlu- 
minosity produced by HBB in more massive models (i.e. with 
M[ > 3.5 - 4.5 M Q , depending on metallicity) is estimated as 
an additional term according to 104 (their equation 32). This 
correction depends on metallicity and it takes into account the 
increase in the core degeneracy during the TP-AGB evolution. 
This means that even when M c remains almost constant in mas- 
sive AGB stars with very deep dredge-up (A « 1), the surface 
luminosity still increases because of the increasing core degen- 
eracy. 

It is important to remark that the M c - L relation is strictly 
valid under the quiescent conditions just before a thermal pulse, 
while significant deviations (i.e. the so-called rapid dip, rapid 
peak and slow dip) take place during the interpulse period due 
to the complex interplay between the occurrence of the flash and 
the consequent structural and thermal readjustment across the 
star. Compared to our previous studies - where a very simpli- 
fied description of the rapid peak and slow dip was adopted - 
a crucial implementation in our synthetic code is that, with the 
aid of the formalism proposed by Wagenhuber & Groenewegen 
(1988) we now follow in detail the evolution of both L(t) and 
Ln(t), which vary significantly during thermal pulse cycles. 

Figures[T]and|2]provide examples of the detailed description 
of the luminosity variations driven by thermal pulses (left pan- 
els) and the evolution of the core mass (right panels) in models 
of different masses and metallicities. These results are in close 
agremeent with full calculations of the TP-AGB phase (see e.g., 
figures 4-9 in Vassiliadis & Wood 1993). We note that, except 
for the final stages of intense mass loss, the most massive mod- 
els (Mi - 4M ) do not show the rapid luminosity peaks because 
of their large envelope mass, and that they are also affected by 



HBB while evolving in luminosity well above the M c - L rela- 
tion. 

It is worth remarking that the growth of the core mass due to 
shell H-burning during the inter-pulses period may be followed 
by a quasi-instantaneous reduction due to convective dredge- 
up at thermal pulses. Such recursive increase-recession of the 
core mass gives origin to the saw-like behaviour displayed by 
all models in the right panels of Figs. Q] and [2] (except for the 
1.0M o , Z = 0.008 case which does not present dredge-up). It 
is already clear from these figures that convective dredge-up is 
predicted to be more efficient in models of larger mass and lower 
metallicity. 

Table Q] presents a few relevant quantities predicted by four 
test synthetic TP-AGB calculations (A-B-C-D models), carried 
out with varying input prescriptions. Model A corresponds to 
our reference case, being obtained with the "best" combination 
of physical assumptions at our disposal. It turns out that the re- 
sults obtained by simply assuming the M c - L relation all over 
the entire interpulse period (model B) are quite different from 
those found with a more realistic description of the pulse cycle 
(model A). In fact, despite the two calculations end up with al- 
most the same final mass, the duration of the whole TP-AGB 
evolution, as well as the M- and C-type phases are significantly 
shorter in model B (by « 30-40 %) than in model A. Model 
C, in which the third dregde-up is suppressed, predicts a size- 
able reduction of the TP-AGB lifetime as well, while the C-type 
phase is completely missed. This simple experiment tells us that, 
for instance, the predicted counts of TP-AGB stars in a galaxy 
via a population-synthesis scheme may be affected by a substan- 
tial error if the underlying TP-AGB tracks are computed using 
a simplified description of the luminosity evolution, i.e. if the 
M c - L relation is assumed anytime. 

2.3. Interpulse period 

The interpulse period T lp marks the pace of thermal pulses dur- 
ing the TP-AGB phase, hence affecting both the evolution over 
the quiescent regime (e.g., core mass and luminosity increase) 
and the events associated to the occurrence of He-flashes (e.g., 
frequency of dredge-up episodes). 

In the present calculations we adopt the formalism sug- 
gested by 104, which is a modified version of Wagenhuber & 
Groenewegen (1988) core mass-interpulse period relation that 
accounts for the increase of Tj p seen in the more massive models 
with deep dredge-up. The 104 (their equation 28) improvement 
consists in adding a positive term, depending on the dredge-up 
efficiency A, to the Wagenhuber & Groenewegen's relation. 

2.4. Effective temperature 

The effective temperature of our TP-AGB models is derived with 
the aid of complete integrations of static envelope models, ex- 
tending from the photosphere down to the core. The adopted in- 
tegration scheme is fully described in Marigo et al. (1996, 1998, 
1999), and Marigo (2002) to whom the reader should refer for 
all details. The fundamental output of the envelope integrations 
is the effective temperature, which affects the properties of the 
TP-AGB evolution, like position in the HR diagram, mass loss 
efficiency, pulsation, colours, etc. 

One crucial aspect to be remarked regards the radiative opac- 
ities at low temperatures (i.e. T < 10 000 K), which are predom- 
inantly due to the molecules (and dust grains for T < 1 500 K) 
that form in the coolest photospheric layers of AGB stars. 
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Fig. 1. Left panels: Evolution of surface luminosity along the TP-AGB phase of two representative models with different initial 
masses with Z = 0.008, according to the analytic formalism by Wagenhuber & Groenewegen (1998). Note the high detail in the 
description of the luminosity variations driven by thermal pulses (i.e. peaks and dips), as well as of the overluminosity caused by 
the occurrence of HBB in the 4.0 M Q model. Right panels: Evolution of the core mass. The saw-like trend reflects the occurrence 
of dredge-up events at thermal pulses. Note that the 4.0 M Q model is characterised by very deep dredge-up for most of its evolution 

Table 1. Predicted properties of synthetic TP-AGB models with Mi = 1.8M , Z = 0.008 under various physical assumptions, 
namely: detailed description of the luminosity evolution (L(t)) or use of the quiescent core mass-luminosity relation (M c - L rel.); 
adoption of molecular opacities for solar-scaled chemical mixtures (ks%) or coupled with the current envelope composition (/c V ar); 
inclusion or not of the third dredge-up. The other table entries correspond to predicted stellar lifetimes (in Myr) covering the whole 
TP-AGB phase (r TP -agbX the M-type phase (tm) with C/O < 1 and log L/ L Q > 3.3, the C-type phase (tc) with C/O > 1 ; final stellar 
mass (Mf) in solar units and surface C/O ratio; effective temperature (in K) at the onset of the superwind regime of mass loss. 
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Importantly, we abandon the unproper choice - still commonly 
used in evolutionary calculations of the AGB phase - of opacity 
tables for solar-scaled chemical composition The present 
TP-AGB models are computed with the adoption of the routine 
developed by Marigo (2002), that calculates the low-temperature 
molecular opacities consistently coupled with the current chem- 
ical composition of the gas (K vax ). In this way we are able to 
account for the drastic opacity changes that take place as soon 



as a model passes from an oxygen-rich (C/O < 1) to a carbon- 
rich (C/O > 1) surface composition as a consequence of the third 
dredge-up, or when the reverse transition possibly occurs due to 
HBB. 

Figures [3] and [4] display the evolution of the effective tem- 
perature and surface C/O ratio during the whole TP-AGB evolu- 
tion of a few selected synthetic models with varying initial mass 
and metallicity. Thermal pulses cause the quasi-periodic com- 
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Fig. 2. The same as in Fig.Q] but for initial metallicity Z = 0.001. 



plex behaviour of T e s that essentially mirrors that of L, shown 
in the corresponding Figs. Q]and[2] As just mentioned, the in- 
troduction of variable molecular opacities is responsable for the 
sudden change in the average slope of the log T e ff (f) curve as 
the photospheric C/O rises above unity in models where the 
third dredge-up take place. In contrast, models without the third 
dredge-up (e.g., M\ — 1 .0 M e , Z = 0.008) are characterised by a 
steady decrease of T e ff for most of their evolution (but for the lat- 
est stages of heavy mass loss). An interesting case is shown by 
the Mi = 4.0 M ,Z = 0.008 model experiencing the sequence 
C/0< 1 — >C/0> 1 — >C/0 < 1 as a consequence of the com- 
petition between the third dredge-up and HBB. Note the kind of 
valley in the T e g curve specularly corresponding to the part of 
the C/O curve above unity. 

As already demonstrated by Marigo (2002) and Marigo et 
al. (2003), the adoption of chemically-variable molecular opaci- 
ties brings a radical improvement in the treatment of carbon-star 
evolutionary models, leading solve several long-lasting discrep- 
ancies between theory and observations. For instance, the red tail 
drawn by field carbon stars in near-infrared colour-colour dia- 
grams of the Magellanic Clouds (e.g., DENIS, 2MASS surveys) 
is reproduced, as well as the low C/O and T e ff values typically 
found in Galactic AGB C-type stars. 

In fact, the reproduction of these basic properties of C stars 
is not possible when using molecular opacities de-coupled from 
the actual surface chemical composition of C-rich models, which 
is still a common choice in published TP-AGB evolutionary cal- 



culations. The sharp dichotomy between the two cases (k^ and 
K V ar) is exemplified in Fig. [5] that compares the predicted be- 
haviour of r eff in a Mi = 1.8, Z = 0.008 TP-AGB model with 
chemically-variable (model A with K val in Table [TJ and with 
fixed, solar-scaled (model D with k^ x in Table[T|i molecular opac- 
ities. It is well evident that the T e s evolution of the model is 
not influenced at all by the behaviour of the C/O ratio, while in 
the K var case the rate of T e $ decrease with time becomes signif- 
icantly faster as the star enters the C-rich domain. As a conse- 
quence, keeping the other input prescriptions the same, the du- 
ration of the TP-AGB phase turns out shorter in the K wm model 
compared to the the Kfi X one. 

Indeed, molecular opacities heavily affect the structure of the 
star's surface layers, as illustrated in Fig. [6] The significant dif- 
ferences between k^ x and K vdl - (top panel) strongly influence the 
temperature (middle panel) and density (bottom panel) stratifica- 
tion in the external regions. In the outermost layers, in particular 
(with mass coordinate 0.998 SS M(r)/M < 1), /e var overcomes 
/Cfix, which determines a sizeable reduction of the effective tem- 
perature, passing from T e ff ~ 328 1 K for to T e ff ~ 2729 K 
for KVar- As further consequence, the A: var model corresponds to a 
more expanded photosphere (i.e. R is forced to increase to allow 
the model radiate the same L ), hence to lower densities (bottom 
panel). 

At this point it is important analysing how much such be- 
haviour is sensitive to the metallicity. Considering that the for- 
mation of molecules in the atmospheres of AGB stars becomes 
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Fig. 3. Evolution of effective temperature and photosperic C/O 
ratio for the same models as in Fig. Q] Note the specular be- 
haviour compared to that of the luminosity, and the remarkable 
average decrease in as soon as C/O increases from below 
to above unity. The reverse trend, i.e. an increase in T e g, occurs 
instead in the 4.0 M Q model due to the subsequent reconversion 
from C/O > 1 to C/O < 1 caused by HBB. 



less efficient at lower Z - because of the higher temperatures 
and the lower abundances of the involved isotopes - it follows 
that the cooling effect on the atmospheres of carbon stars due to 
molecular opacities should be less and less evident at decreas- 
ing metallicity. To give an example, by comparing Fig. [3] and 
Fig. [4] one already notices that the sudden change in the mean 
slope of the effective temperature curve, at the transition to the 
C-star domain, is much more pronounced in TP-AGB models 
with Z = 0.008 compared to those with Z = 0.001. This point 
is better illustrated with the help of Fig. [7] that shows how the 
effective temperature of a selected model (with given total mass, 
core mass, and luminosity) is expected to behave at increasing 
C/O ratio over a wide range of initial metallicities. As already 
mentioned, the sudden and extended excursion toward lower ef- 
fective temperature as soon as C/O^ 1, characterising models 
of higher Z, is progressively reduced at decreasing metallicity, 
eventually disappearing for Z < 0.0004. 

2.5. The third dredge-up 

The treatment of the third dredge-up is essentially the same one 
as in Marigo et al. (1999) to which the reader should refer for all 



Fig. 4. The same as in Fig. [3] but for initial metallicity Z = 0.001 . 



details. In practice we need to specify three quantities, namely: 
(i) the minimum temperature at the base of the convective en- 
velope at the stage of the post-flash luminosity maximum, r* ed , 
or equivalently the minimum core mass required for the occur- 
rence of the mixing event, M™ n ; (ii) the efficiency of the third 
dredge-up A; and (iii) the chemical composition of the inter- 
shell. Compared to our previous models (e.g., Marigo & Girardi 
2001 and references therein), in this work we improve the de- 
scription of the third dredge-up, thanks to the recent results of 
detailed calculations of the AGB phase, namely K02 and 104. 
The main points are summarised as follows. 

2.5.1. Occurrence of dredge-up 

In Marigo et al. (1999) the adopted criterion assumes that t£ kA , 
derived from numerical integrations of static envelope models, 
is constant regardless of stellar mass and metallicity, following 
the indications from full AGB calculations as reported by Wood 
(1981). The calibration of the parameter from the fits to the ob- 
served CSLFs in the Magellanic Clouds yields log T* ed = 6.4. 

In this study we explore an alternative approach to log 7^ red 
by introducing the quantity M™ n , that is the minimum core mass 
required for dredge-up to take place. This has been a commonly 
used parameter in synthetic AGB models over the years, often 
assumed to be the same (~ 0.58 M ) for any stellar mass and 
metallicity (Groenewegen & de Jong 1993; Marigo et al. 1996; 
Mouhcine & Lancon 2002). Only very recently extended grids of 
complete AGB models have become available (e.g., K02; 104), 
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Fig. 5. Predicted behaviour of the effective temperature (top 
panel) and photospheric C/O ratio (bottom panel) over the en- 
tire TP-AGB evolution of a model with Mi = 1.8M , Z = 0.008, 
adopting either fixed solar-scaled (model A in Table[TJ gray line) 
or variable (model D in Table [T] black line) molecular opaci- 
ties. The two r e ff curves suddenly separate as soon as the mod- 
els becomes carbon-rich as a consequence of the third dredge- 
up. Arrows indicate the quiescent stage just preceeding the 24 th 
thermal pulse in both models, whose photospheric structures are 
shown in Fig. [6] 

allowing to relax the crude simplification of constant M™ n in 
place of a more physically sound dependence on stellar mass 
and metallicity. 

In the present study we adopt the relations M™ n (M, Z), pro- 
posed by K02, based on grids of complete AGB calculations 
covering initial masses M\ — 1 - 6M and metallicities Z 
0.02. 0.008, 0.004. A partial modification of their analytical fits 
turns out necessary for the following reasons. 

Systematic shift. The reproduction of the observed CSLFs in 
the Magellanic Clouds requires that M™ n is lower than predicted 
(Sect. l4~Tl see also K02 and 104). Our synthetic calculations in- 
dicate that M™ n - as given by equation 10 in K02 - must be 
decreased by an amount AM™ n that depends on metallicity: 
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Fig. 6. Structure of the surface layers of a model with M = 
1.776M , M c = O.632M , logL = 3.926L and C/0= 2.560, 
as a function of the mass coordinate M{r)/M, where M(r) is 
the mass within a sphere of radius r. From top to bottom pan- 
els we compare the behaviour of the Rosseland mean opacity 
(cm 2 gr~'), the temperature, and the density (in g cm -3 ) for 
(solid line) and /c val (dashed line) cases. 



The negative correction, -AM™ 11 , becomes more important at 
decreasing metallicity, being e.g., (-0.01, -0.05, -0.08) M for 
Z = (0.008, 0.004, 0.001) respectively. This implies that for Z < 
0.001 the onset of the third dredge-up in stellar models of any 
mass would take place since the first thermal pulse. 

In practice, we limit this metallicity extrapolation imposing 
that no dredge-up can occur for initial masses M < 1 M , in 
order to comply with the observed lack of intrinsic C-stars (not 
related to binary evolution like the CH stars or dwarf carbon 
stars) in Galactic globular clusters. 

Correction for overshooting As far the evolution prior to the 
AGB is concerned, K02 do not assume overshoot from convec- 
tive cores during H-burning, while a moderate convective over- 
shoot is included in Girardi et al. (2000) sets of tracks, from 
which we extract the initial conditions of our AGB models. 

One important effect due to convective overshoot is that of 
lowering the critical mass Me e F, that is maximum initial stel- 
lar mass required to ignite core He-burning in degenerate condi- 
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Fig. 7. Metallicity dependence of the relationship between the 
effective temperature and the photosphetic C/O ratio. The curves 
are obtained by means of envelope integrations assuming the 
same set of stellar input parameters (i.e. total mass, core mass, 
luminosity), while varying the initial metallicity, as indicated. 
Starting from solar C/O, the carbon abundance is progressively 
increased to mimic the effect of the third dredge-up. We notice 
that the sudden photosperic cooling caused by molecular opaci- 
ties as soon as C/O > 1 is quite pronounced at higher metallici- 
ties, while it almost disappears for Z < 0.0004. 

tions. At given metallicity this limit mass is well defined by the 
minimum in the relation between the core mass at the first ther- 
mal pulse, Ml, and the initial stellar mass M; (see e.g., Marigo et 
al. 1999; K02). A similar behaviour is shared by the correspond- 
ing M™ 11 vs M[ relation (see e.g., K02), so that the minimum of 
M™ n corresponds to M He F- 

It follows that the inclusion of convective overshoot in stellar 
evolutionary calculations causes, among other effects, a decrease 
of M^f, hence a shift of the minimum M™ n (Mi) towards lower 
masses. In other words, the inclusion of convective overshooting 
can be approximately described by considering that a star with 
initial mass M\ behaves like a star of higher mass f ow M\ (with 
/ov >1) without overshooting. 

From a comparison between K02 and Girardi et al. (2000) 
models the enhancement factor is estimated f m ~ 1.2. Hence, in 
order to mimic the effect of convective overshoot on the relation 
M™ n (Mj) we simply adopt the K02 formalism as a function of 
the adjusted mass variable M = f ov M: 



Mf n (M) = [jlC°(Af)] 



K02 



(2) 



The resulting M™ n at varying initial stellar mass are shown in 
Fig. E]for Z = 0.019, 0.008, 0.004. For metallicities outside the 
range covered by K02 models, we linearly extrapolate in Z the 
coefficients of the fitting polynomial functions of K02. It should 
be recalled that for models with initial mass M\ <: 3M G the core 
mass predicted by Eq. (f2| is lower than the core mass at the 
first thermal pulse M C; itp, so that in these cases we set Mf m = 

M C; iTP- 
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Fig. 8. Minimum core mass required for the occurrence of the 
third dredge-up as a function of the stellar mass at the first ther- 
mal pulse, for three choice of the initial chemical composition. 
Note that in all curves the minimum corresponds to the critical 
mass M He . See text for more explanation. 



2.5.2. Dredge-up efficiency 

We briefly recall that the efficiency of each dredge-up event is 
usually expressed via the quantity A - AMd re d/AM c , i.e. the 
mass brought up to the surface at a thermal pulse, AMd re d, nor- 
malised to the core mass increment over the preceeding inter- 
pulse period, AM C . 

The usual choice in previous studies has been adopting the 
same A for any value of the stellar mass. The past calibration of 
this parameter through the reproduction of the observed CSLFs 
in the Magellanic Clouds (see e.g., Marigo et al. 1999) yields 
A = 0.5 for Z = 0.008 and A = 0.65 for Z = 0.004. 

In the present study we relax the assumption of constant A, 
accounting for both mass and metallicity dependences found on 
the base of K02 and 104 complete AGB models. In summary, 
once M c has grown above M™ n , the efficiency A progressively 
increases towards an asymptotic value A max , with an e-folding 
interval of N r thermal pulses. This behaviour is conveniently de- 
scribed by the relation 



M.N) = imax [1 - SX P (N/N r )] 



(3) 



where N is the progressive pulse number (N — as long as 
M c < Mf n ). 

As to /l max we use the analytical fits by K02. We notice that at 
given metallicity the asymptotic quantity increases with the stel- 
lar mass, typically reaching values 0.8-1.0 for Mj > 3 M . The 
metallicity dependence is less strong, such that at given stellar 
mass A max tends to be somewhat higher at decreasing Z. 

In order to have significantly larger efficiencies at metallici- 
ties Z < 0.004 - that is below the low end of the metallicity range 
covered by K02 models - we modify the original A m . dx (M, Z) re- 
lation in K02 by applying a positive shift to the mass variable 



^max(M) = A mm (M + AM A ) 



(4) 
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Table 2. Coefficients for Equation (O 







Z 




0.019 


0.008 


0.004 




4.110 


2.785 


2.555 


a 2 


42.612 


10.625 


809.426 


«3 


5.834E-03 


4.806E-02 


3.348E-05 


"4 


-2.113 


-0.908 


-4.417 


Cls 


2.014 


0.806 


1.751 


a h 


-9.116 


-6.708 


-70.624 


a~i 


3.830 


3.462 


3.516 



where 



AM A = 200 (0.004 - Z) M e 



(5) 



This correction is consistent with the theoretical expectation 
that the third dredge-up is more efficient at lower metallicities 
(see section 3.4 in 104 for a similar adjustment), and it is also 
required to reproduce the faint wing of the CSLF in the SMC, 
which is populated mainly by lower-mass stars with metallicities 
down to * 0.001 (see Sect.©. 

As discussed by K02 there is no simple law able to fit N T as 
a function of M and Z. A rather satisfactory reproduction of the 
results for N T presented by K02 (in their table 5) is obtained with 
the function 



iV r = d\ + «2 [«3 - exp(fl4Mixp)] 

X [g^ + 06 exp(fl7 - Mitp) 2 ] . 



(6) 



The coefficients a,- , i = 1 , . . . , 7 are given in Table [2] for Z = 
0.019, 0.008, and 0.004. Interpolation and extrapolation in logZ 
are adopted for other metallicities. 

2.5.3. Chemical composition of the inter-shell 

The abundances of 4 He, 12 C, and le O in the intershell region 
are derived as a function of the core mass growth AM C = M c - 
M Cj itp on the base of nucleosynthesis calculations by Boothroyd 
& Sackmann (1988). A fit to their figure 9 gives: 



X( 4 He) 
X( 12 C) 
X( 16 C) 



0.95 + 400 (AM C ) - 20.0 AM C 
0.03 - 352 (AM C ) 2 + 17.6 AM C 



(7) 



-32(AM C ) Z + 1.6 AM C 

for AM C < 0.025 M Q , and 

X( 4 He) = 0.70 + 0.65 (AM C - 0.025) 
X( l2 C) = 0.25 - 0.65 (AM C - 0.025) 
X( 16 0) = 0.02 - 0.065 (AM C - 0.025) 



(8) 



for AM C > 0.025 M Q . These relations are displayed in Fig. [9] 

Looking at Fig. [9] we notice that the fractional mass of pri- 
mary I2 C quickly increases during the first pulses up to a max- 
imum value of 0.25 for AM C = 0.025 M Q , and then it starts to 
decrease slowly. The same behaviour characterises the evolution 
of the I6 abundance (with a maximum value of 0.02), while 
a mirror-like trend is followed by 4 He that reaches a minimum 
value of 0.70. 



2.6. Hot-bottom burning nucleosynthesis 

HBB takes place in the most massive TP-AGB models, with 
(Mi > 3.5 -4.5M ) depending on metallicity, and it corresponds 
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Fig. 9. Chemical composition ( 4 He, 12 C, and 16 abundances in 
mass fraction) of the dredged-up material as a function of the 
core mass increment during the TP-AGB phase, as predicted by 
Eqs. and ©. 
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Fig. 10. Evolution of the temperature at the base of the convec- 
tive envelope during the whole TP-AGB phase of the (4M , Z = 
0.008) model that experiences both the third dredge-up and 
HBB. The inset zooms the temporal behaviour of T\, ce during 
few selected pulse cycles. 



to H-burning, mainly via the CNO cycle, in the depeest layers 
of their convective envelopes. For the present work we adopt the 
same treatment of this process as in Marigo et al. (1999) to which 
we refer for more details. 

From an observational point of view, HBB does not only 
make models more luminous than expected by the core mass- 
luminosity relation (see Sect. 12.21 ) but, more importantly, it pro- 
duces also notable chemical changes in the surface chemical 
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composition. One of the most significant and well-known effects 
played by HBB nucleosynthesis is that of converting C into N, 
hence delaying or even preventing the formation of C stars. 

However, it should be noticed that, under particular circum- 
stances, HBB might instead favour the formation of C stars 
thanks to the destruction of O in favour of N. This applies to 
low-metallicity massive models in which the temperature at the 
base of the envelope is high enough to activate the ON cycle. 
A few exemplifying cases of models experiencing both the third 
dredge-up and HBB are discussed in Sect. [33] and illustrated in 

Fig.EU 

It is worth emphasizing that, at each time step during the in- 
terpulse periods, the elemental abundances in the envelope of 
TP-AGB models with HBB are determined by solving a net- 
work of nuclear reactions (the p-p chains and the CNO tri- 
cycle), according to the current temperature and density strat- 
ifications given by numerical integrations of a complete enve- 
lope model (see Marigo et al. 1998 for a complete description 
of the method). This allows a consistent coupling of the HBB 
nucleosynthesis with the structural evolution of the envelope, a 
technical characteristic that places our synthetic TP-AGB model 
closer to full TP-AGB models. In fact, as already shown by 104 
a parameterised treatment of HBB nucleosynthesis in synthetic 
TP-AGB calculations turns out quite complex and sometimes 
troublesome for the relatively large number of the parameters 
involved (« 7), their intrisic degeneracy, and the non-univocity 
of the calibration. 

As an example let us consider one of these quantities that 
is directly related to the efficiency of HBB, i.e. the temperature 
at the base of the convective, T\, ce . 104 (see their section 3.3.4) 
propose a parametric formalism for 7j, ce , which in turn requires 
additional parameters to be calibrated on the base of full TP- 
AGB calculations. 

On the other side, by means of envelope integrations car- 
ried out with our synthetic TP-AGB code, we are able to natu- 
rally predict the temporal evolution of 7b ce during the whole TP- 
AGB phase, as displayed in Fig. [10] for the (Mj = 4.OM , Z = 
0.008) model. Note the complex behaviour and wide excursion 
in 7b ce due to thermal pulses, which are zoomed in the inset. 
Considering only the quiescent pre-flash stages, that ideally de- 
scribe the upper envelope of the 7j,ce curve, we see that the tem- 
perature steeply rises during the first pulse cycles, then it flattens 
out toward a maximum value, and eventually it drastically drops 
as soon as the envelope mass is significantly reduced by mass 
loss. 

2.7. Pulsation modes and periods 

Models for pulsating AGB stars (see e.g., Fox & Wood 1982) in- 
dicate that the dominant growth rate (GR) among different pul- 
sation modes is correlated with luminosity, i.e. modes of higher 
order have the highest GR at lower luminosities. This means 
that if we consider, for instance, the fundamental (FM) and first 
overtone (FOM) modes with growth rates GRo and GR], respec- 
tively, we expect GRi >GRo at lower luminosities while the sit- 
uation is reversed at higher luminosities. 

Hence making the reasonable assumption that the active pul- 
sation mode among the possible modes is that one corresponding 
to the highest GR, we derive a picture of pulsation on the AGB: 
as a star climbs in luminosity it becomes unstable through a se- 
quence of pulsation modes, first exciting those of higher order 
(i.e. third — > second — > first overtone) and eventually switching 
to the fundamental mode (see Lattanzio & Wood 2003). 



Due to the lack of sufficient theoretical data on higher order 
pulsation modes, in this study we limit to consider the first over- 
tone mode (FOM) and the fundamental mode (FM). We assume 
that all stellar models start their TP-AGB evolution as FOM pul- 
sators. The transition from the first overtone to the FM pulsation 
may take place later, provided a minimum luminosity, Li_o, is 
reached. This critical luminosity is defined by the stage at which 
the growth rates per unit time, GR, for both modes become equal 
(GRi =GR fl 

To obtain Lj_o we use the results of the linear, nonadiabatic 
pulsation models for Mira variables calculated by Ostlie & Cox 
(1986). A good fitting relation (based on their table 2) providing 
the critical luminosity as a function of stellar mass and effective 
temperature is 

logLi_o/£o = -14.516 + 2.277 log M/M (9) 
+5.046 log r eff . 

As we notice from Fig.[TT]the switch from the FOM to the FM 
is expected to occur at higher luminosity for larger stellar mass 
and effective temperature. Since Ostlie & Cox (1986) present 
models for solar composition only - with (T, Z) = (0.28, 0.02) -, 
we include a metallicity dependence to Eq. (0 as quantified by 
Bedijn (1988) on the base of pulsation models for long-period 
variables calculated by Fox & Wood (1982) for metallicities Z = 
0.001,0.01,0.02. Following Bedijn (1988, see his equation 22) 
the transition radius, Ri-o, from the FOM to the FM scales as 
(Z/0.02)- 042 , hence we get 

logL^o/io = -14.516 + 2.277 log M/M (10) 
+5.046 log r eff - 0.084 log(Z/0.02) , 

considering that L;_o K ^i-o^Tff ^ or gi ven effective temperature. 

The pulsation periods are derived from theoretical (P, M, R)- 
relations for long-period variables calculated by Fox & Wood 
(1982). Specifically, FOM periods (Pi) are derived through the 
analytical fit of the pulsation constant (Wood et al. 1983): 

g/days = (P I /days)(M/M ) 1/2 (/?//? )- 3/2 (11) 
= 0.038 

(if M > 3M ; or M = 2M and Pi < 300 days; 
or M = 1M and Pi < 200 days) 

= 0.038 + 2.5 x 10~ 5 (Pi - 300) 
(if M = 2M and P 1 > 300 days) 

= 0.038+ 4.5 xl0~ 5 (Pi- 150) 
(if M = IM Q and P { > 150 days) 

= 0.038 + 5.5 x 10~ 5 (Pi - 100) 
(if M = O.7M and P { > 100 days) 

whereas FM periods, Pq, are calculated with 

log(/Vdays) = -2.07 + 1.94 log(R/R e ) (12) 
-0.9 log(M/M ) (if M < 1.5M ) 
= -2.59 + 2.2 log(R/R Q ) 

-0.83 log(M/M ) - 0.08 log(Z/10~ 3 ) 
+0.25(F - 0.3) (if M > 2.5M ) 

where Z and Y denote metallicity and helium content (in mass 
fraction), respectively. These latter formulae are taken from 

2 This latter criterion is the same one as proposed by Bedijn (1988) 
in his synthetic AGB model. 
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Fig. 11. Critical luminosity Li_o defined by the equality of the 
growth rates for FO and FM pulsation. Symbols show the predic- 
tions of nonadiabatic pulsation models by Ostlie & Cox (1986) 
for different values of stellar mass and effective temperature. 
Lines (a-e) are plotted according to the fitting Eq. (TTOt . each of 
them corresponding to a given T e ff at varying stellar mass. See 
text for more details. 



Groenewegen & de Jong (1994; for M < 1.5M ) and Fox & 
Wood (1982; for M > 2.5M ). For 1.5 < (M/M ) < 2.5, logP 
is linearly interpolated using log M as the independent variable. 

2.7.1 . A cautionary remark: how to treat pulsation in C-rich 
models 

At present all published pulsation models for long-period vari- 
ables (e.g., Fox & Wood 1982; Wood 1990; Ya'Ari & Tuchman 
1996) are strictly valid only for oxygen-rich chemical compo- 
sitions, with photospheric C/0< 1, while models suitably con- 
structed for pulsating C stars - and using proper C-rich opacities 
- are still missing. This means that the extension of the (P, M, R)- 
relations available in the literature, is not straightforward, and 
possibly wrong. Since one of the main goals of the present study 
is to provide a more reasonable description of the AGB evolution 
as a function of the stellar chemical-type, we face the problem 
to assign periods to C-rich models. Given the lack of theory, the 
only remaining approach comes then from the empirical ground. 

For instance, assuming that the inertial radius R, which de- 
termines the pulsation period, coincides with the photospheric 
radius - defined by the black-body relation L = R 2 T^ (in so- 
lar units) -, and applying Eqs. ( [TTIfT2l to compute periods of 
C-rich models, then one would expect that, at a given luminos- 
ity, C stars should have larger R owing to their lower T e ff , hence 
longer periods compared to M stars. 

As already noted long ago by Wood et al. (1983) this predic- 
tion is contradicted by observations since both groups of stars 
appear to populate the same sequences in the (Mboi - P) dia- 
gram, as shown for instance in figure 6 of Wood et al. (1983; see 
also Hughes & Wood 1990). This feature has been confirmed by 
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Fig. 12. Evolution of pulsation period along the TP-AGB phase 
of a few selected models with initial metallicity Z = 0.008. First- 
overtone and fundamental periods, P - P\ and P = Pq, are 
drawn with thick and thin solid lines, respectively. 



more recent and extended surveys of long-period variables (e.g., 
Groenewegen 2004; Fraser et al. 2005). 

As suggested by Wood et al. (1983) the immediate implica- 
tion is that, differently from the photospheric radius, the inertial 
radius of pulsating AGB stars should be largely insensitive to the 
photospheric C/O ratio. In other words, while the photospheric 
radius of C stars is expected to be much more extended than that 
of M stars of the same luminosity due to an opacity effect, the 
inertial radii of both classes of stars should be quite similar. 

For all these reasons, in this work we simply decide to cal- 
culate periods as follows (all quantities are in solar units) 

R = L° 5 r e 2 ff for M-type models with C/O < 1 (13) 

that is the inertial radius coincides with the photospheric radius, 

R = L° 5 r e 2 ff M for C-type models with C/O > 1 (14) 

that is the inertial radius is calculated by assigning the effective 
temperature T e ^ _m that would compete to an M-type model with 
the same luminosity. In these case is derived from enve- 

lope integrations always assuming a fictious oxygen-rich com- 
position (i.e. assuming the solar ratio, C/0= 0.48). 

The predicted period evolution during the whole TP-AGB 
phase of a few selected models is shown in Figs. [12] and [13] 
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Fig. 14. Mass evolution of O-rich TP-AGB models with Z = 
0.019 as a function of luminosity, with mass-loss rates cal- 
culated on the base of Bowen & Willson (1991) dynamical 
atmospheres for Miras including dust. The five parallel solid 
lines define approximate loci of constant M/M (in units of 
yr -1 ), corresponding to evolutionary slopes dlogM/dlogL = 
-0.0018,-0.018,-0.18,-1.8,-18 (from left to right). The 
dashed line is the "cliff" line, along which dlogM/dlogL = 1 
and M/M = 5.67 x 10~ 7 yr -1 . 



Fig. 13. The same as in Fig. [12] but for metallicity Z = 0.001. 



2.8. Mass loss rates 

Mass loss rates, M, during the TP-AGB evolution are calculated 
in a different fashion depending on whether a stellar model is 
oxygen-rich or carbon-rich. 



The variations of luminosity and photospheric radius produced 
by thermal pulses reflect in the quasi-periodic behaviour of P, 
in substantial agreement with the predictions of detailed AGB 
models and their pulsation properties (e.g., Vassiliadis & Wood 
1993; Wagenhuber & Tuchman 1996). Moreoever, in all models 
it is well evident the switch from FOM to FM pulsation regime 
and viceversa, corresponding to a mean change in period by a 
factor of ~ 2. According to the luminosity criterion given by 
Eq. ( [Tol l the mode switching would not correspond to a single 
event during the TP-AGB evolution, but it could actually take 
place several times, driven by the luminosity variations produced 
by thermal pulses. As displayed in Figs. [12] and [13] there might 
be cases in which the mode transition may occur during the pulse 
cycle, with the FOM pulsation charactering the low-luminosity 
dips, and the FM pulsation being recovered at increasing lumi- 
nosities during the quiescent stages. It should be also noticed 
that the r e ff-dependence of Li_o in Eq. ( fTOb favours the switch- 
ing to the FM pulsation at decreasing effective temperature. This 
explains the occurrence of a short FM stage experienced by the 
(Mi = 4.OM , Z = 0.008) model at evolutionary ages t « 5 x 10 5 
yr, which just corresponds to the transitory fulfillment of the 
condition C/O > 1 (see Fig. 0), thanks to the effect of the third 
dredge-up prevailing over that of HBB. 



2.8.1. O-rich models 

In the C/O < 1 case we use the results of dynamical atmospheres 
including dust for long-period variables calculated by Bowen & 
Willson (1991). The grid of model atmospheres is suitably con- 
structed for oxygen-rich Mira-type stars with solar metallicity 
and it covers the ranges of initial masses M\ = 0.7 - 2.4 M and 
periods P — 150 - 800 days. Pulsation is assumed to occur in the 
FM. 

By suitably combining their dynamical calculations and with 
a few basic formula^] to describe the evolution of stellar parame- 
ters on the AGB, Bowen & Willson (1991) derive a relationship 
between the mass loss rate M and the slope of a stellar AGB 
evolutionary track in the log M-log L diagram: 



^=5.65xKT 7 ^I^ 
M dlogL 



[yr -1 ] 



(15) 



where M and L are in solar units, and M is given in M Q yr -1 . 

To summarise, the adopted procedure to specify M for an 
oxygen-rich AGB model is as follows. Given the current values 



3 These refer to the period-mass-radius relation for FM pulsation, the 
radius-luminosity-mass relation, and the core mass-luminosity relation. 
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of stellar mass M, photospheric radius R, and pulsation period 
P (obtained according to Sect. 12.71 ), we derive the corresponding 
luminosity used by Bowen & Willson (1991, their equation 3) in 
their computation of model atmospheres. Then, for any specified 
combination of (L, M), the evolutionary slope dlogM/dlogL 
is obtained by interpolation between grid points defined by the 
lines of constant M/M shown in figure 2 of Bowen & Willson 
(1991). Finally the mass-loss rate follows from Eq. (15[ . 

Figure [14] shows a grid of TP-AGB evolutionary tracks with 
C/O < 1 in the log M-log L diagram. We plot only models at the 
quiescent stages of pre-fiash luminosity maximum. Along each 
track the stellar mass keeps constant until an abrupt change of 
slope takes place. This marks the onset of the superwind mass 
loss, which determines the quick ejection of the whole envelope 
and the termination of the AGB phase. The sequence of these 
sharp corners at varying M and L form the so called "cliff" line, 
as designated by Willson (2000). It is defined by the condition 
d log M/dt ~ d log L/dt and it corresponds to a critical mass loss 
rate M crit ~ 5.67 x 10~ 7 x M in units of M yr _1 . Once reached 
the "cliff" point the star enters the superwind phase. 

As anticipated by Bowen & Willson (1991) and thoroughly 
discussed by Willson (2000), at decreasing metallicity the mass- 
loss rate of an oxygen-rich AGB star is affected by two fac- 
tors, namely: i) less efficient dust production due to the lower 
abundances of involved elemental species, and ii) more compact 
structure, i.e. smaller photospheric radius at given luminosity. It 
turns out that the latter effect is the largest, producing a shift of 
the cliff line to higher L for lower Z (see figure 9 in Willson 
2000). Following Bowen & Willson (1991) for Z < 0.1Z o the 
onset of the supewind takes place at a critical luminosity, ap- 
proximately given by 



^sw(Z) - Lsw(Z ) + ALsw, 
where 

ALsw = 0.12-0.13 log(Z/Z G ) 



(16) 



(17) 



A few comments should be made at this point. First, we 
should keep in mind that the results obtained by Bowen & 
Willson (1991) and illutrated in the key (logM vs logL) dia- 
gram (their figure 2) are not only dependent on the details of 
their model for dynamical atmospheres, but they also rely on the 
adopted relations to describe the underlying AGB evolution for 
various choices of the stellar mass. This fact introduces some in- 
consistency in our calculations as we are dealing with different 
synthetic TP-AGB models (e.g., different core mass-luminosity 
relation). Nevertheless the adopted approach is the best possible 
with the available theoretical information. Moreover, it produces 
results that considerably improve the comparison between pre- 
dictions and observations (see Sect. 13.51 ). 

Second, we recall that Bowen & Willson (1991) calculations 
are strictly valid for FM pulsation. There is no extended anal- 
ysis of the same kind devoted to pulsation in modes of higher 
order. A limited set of dynamical models for FOM pulsation in- 
cluding dust are presented by Bowen (1988). By comparing the 
mass-loss rates calculated by Bowen (1988) for FM (table 6) 
and FOM models (table 8) with the same given combination of 
stellar mass, radius and effective temperature^ we obtain the be- 
haviour shown in Fig. [17] 



4 For assigned fundamental-mode periods P = 350, 500, 700 days, 
and stellar masses, M = 0.8, 1.0, 1.2, 1.6, 2.0 M Q , the corresponding 
radii and first-overtone mode periods are derived from the period- 
mass-radius relationships by Ostlie & Cox (1986), while it is assumed 
r cf f = 3 000 K in all cases. 



io- 



io- 



10" 



10- 4 



fc; 10-° 



1.0 M Z='0.008 



5xl0 5 



10 6 



1.5xl0 6 



10" 



10" 



r'1.8 M Q k = 0.008 



10 6 



2x10 s 3xl0 6 4xl0 6 




2xl0 5 4x10 s 6xl0 5 8xl0 5 10 6 

t (yr) 

Fig. 15. Evolution of the mass-loss rate along the TP-AGB phase 
of a few selected models. 

Table 3. Coefficients for Equation ( fTSl 
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3.9290 
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-7.2055 


-6.0280 


C3 
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3.4831 


25.6437 



We notice that the trend is not monothonic, as M for P = 
P\ (hereafter M\) may be lower or larger than M for P — Po 
(hereafter Mo) depending on stellar mass. Moreover, the peak in 
the Mi /Mo ratio appears to become higher and to shift towards 
lower masses at decreasing Pq. 

Polynomial fits of M\/Mq for three values of Po are derived 
assuming the rational form 



Mi bi + biM + b^M 2 
M c i + ciM + c 3 M 2 



(18) 



The coefficients a, and bi are given in Table [3] 

In practice, in order to estimate the mass-loss rate of a TP- 
AGB star pulsating in the FOM mode we proceed as follows. For 
given stellar mass, radius and luminosity we compute the period 
Pq via eq. (TT21 and the corresponding mass-loss rate Mo accord- 
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Fig. 16. The same as in Fig. [15] but for metallicity Z = 0.001. 

ing to Sect. 12.81 Then we apply Eq. ( fl~8l to derive M\ (linear 
interpolation or extrapolation in log P are adopted between grid 
points at given M). 

Being aware that this recipe based on the few FOM models 
calculated by Bowen (1989) does not allow to derive a com- 
plete picture of mass loss under FOM pulsation, nevertheless 
it represents a first attempt to address the question. Moreover, 
this is likely a more realistic approach than simply scaling the 
mass-loss rate for FM pulsators for P = P\ » Pq/2.2 (cf. 
Fig. |23l l, which would always lead to very low mass loss-rates. 
On the contrary, observations indicate that low -period AGB vari- 
ables, like the semiregulars, are characterised by mass-loss rates 
comparable with those of Miras, up to few 1O~ 7 M yr~' (e.g., 
Olofsson et al. 2002). If these stars are assigned a pulsation mode 
different from the FM, likely the FOM, as currently suggested 
(see Lattanzio & Wood 2003), then this would imply that i) ei- 
ther the M(P) relation depends on the pulsation mode, or ii) M 
is independent of P and the observed correlation between M and 
P just hides the true dependence on other basic stellar parame- 
ters, like mass, luminosity and effective temperature. It is clear 
that detailed analyses of dynamical FOM pulsating atmospheres 
should be of paramount importance to investigate this point. 

2.8.2. Carbon-rich models 

In the C/O > 1 case we exploit the results of pulsating dust- 
driven wind models for carbon-rich chemistry, developed by the 
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Fig. 17. Predicted ratio between the mass loss rates M\ and Mo, 
under the assumption of FOM and FM pulsation respectively, 
as a function of the stellar mass according to the dynamical at- 
mosphere models calculated by Bowen (1988). Three groups of 
models, characterised by the same Pq, are shown together with 
the corresponding polynomial fits. 



Berlin group (Fleisher et al. 1992; Sedlmayr & Winters 1997; 
Winters et al. 2000; Wachter et al. 2002; Winters et al. 2003 
and references therein). Extended grids of models are suitably 
constructed for carbon-rich AGB stars, yielding the mass-loss 
rate as a function of main stellar parameters, such as: effective 
temperature (2200-3000 K), luminosity (3000-24000 L @ ), mass 
(0.8-1.8 M e ), photospheric C/O ratio (1.20-1.80), and pulsation 
period (100-1000 days). According to the detailed investigation 
by Winters et al. (2000, 2003) one can schematically distinguish 
two dynamical regimes for the stellar outflows. 

The first regime refers to the so-called B-type models, and 
it should apply to low-luminosity, short-period, high-effective 
temperature red giants, for which the dominant mass-loss driv- 
ing mechanism is stellar pulsation, while radiation pressure on 
dust grains plays a reduced role. In fact for this class of models 
radiation pressure alone is not able to produce an outflow, such 
condition being expressed by a < lQ During this regime tenu- 
ous winds are generated, characterised by low terminal veloci- 
ties (v'oo < 5 km s _1 ), and low mass-loss rates, never exceeding a 
critical value of M clit « 3 x 10~ 7 M yr -1 . 

The subsequent development of a stable dust-driven wind 
from a carbon star (with M > M cl -i t ) requires a minimum lu- 
minosity Lsw, such that the outward acceleration provided by 
radiation pressure on dust exceeds the inward gravitational pull 
(a > 1). According to Winters et al. (2000) models this critical 
luminosity depends mainly on stellar mass and effective temper- 
ature. As soon as L > Lsw, a carbon star is expected to enter 
the dust-driven-wind regime, described by the so-called A-type 
models, with mass-loss rates in the typical range 3 x 10~ 7 % M % 



5 The standard notation of dynamical wind models defines a = 
a md /a grav , the ratio of the radiative acceleration over the gravitational 
one. 
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Fig. 18. Critical luminosity for the onset of the super-wind in C- 
rich models as a function of the initial stellar mass, according 
to Eq. ([Till (dashed line) for P = 400 days, and Eq. @ (solid 
line) for three choices of P, as indicated. In all cases the effective 
temperature is set equal to 2 800 K. 

1O _4 M yr , and terminal velocities in the range 5 ^ v«> S= 40 
km s _1 . 

The picture just depicted leads us to adopt the following 
recipe to calculate the mass-loss rate for carbon-rich models 
evolving on the TP-AGB. First we need to specify the critical 
luminosity Lsw that separates the two wind regimes. Basing on 
wind models calculated by Schroder et al. (1999) - with tran- 
sition properties between the B and A classes - Schroder et al. 
(2003) propose an approximate relation, giving Lsw as a func- 
tion of stellar mass M and effective temperature T e g : 

log(L sw /L ) = 3.8 + 4(log L eff - 3.45) + log(M/M ) (19) 

The above relation is derived from a small subset of wind mod- 
els for various choices of stellar mass and effective temperature, 
while adopting the same period, P = 400 days, and carbon-to 
oxygen ratio, C/0= 1.3. 

In the present work an attempt is made to refine this relation 
by including also a dependence on the pulsation period P. To 
this aim we use a sample of pulsating wind models computed by 
Winters et al. (2000; their tables 6-9) at varying input parame- 
ters, namely: stellar mass M, effective temperature L e ff, pulsa- 
tion periods P (in days), carbon-to-oxygen ratio C/O, and piston 
velocity Av p . 

From these data we derive the stellar luminosity at which 
the time-averaged (over a few pulsation periods) a « 1, as a 
function of M, L e ff, and P (in days), while we consider fixed both 
the carbon-to-oxygen ratio (C/0= 1 .2), and the piston velocity 
(Av p = 5 km s ). The resulting fitting relation reads 

log(L sw /L ) = 3.796 + 6.614 log(J e ff/2800) (20) 
+ 1 .096 log(M/M ) - 1.728 log(P/400) . 

For the sake of comparison Fig. [I8]displays Lsw as a func- 
tion of stellar mass, as predicted for T e g = 2800 K by our pro- 
posed fit (Eq. l20l for three choices of the pulsation period, and 
by Eq. JT9] l without period-dependence suggested by Schroder 
et al. (2003). We note that our proposed formula recovers the 

6 In the piston-approximation, commonly used in pulsating wind 
models, the piston amplitude Ai' p denotes the velocity amplitude of the 
innermost grid point, assumed to vary sinusoidally over the pulsation 
period P. 



Schroder et al. (2003) relation well for P = 400 days, while pre- 
dicting a significant sensitiveness of Lsw to period changes, i.e. 
the onset of the dust-superwind should occur at lower luminosi- 
ties for longer periods. 

As long as L < Lsw the mass-loss rate is calculated as a 
function of the near infrared J - K color 

logtM/tMoyr- 1 )] = - ~ 6 ° -4.0 (21) 

J — A + U.l 

according to the semi-empirical calibration by Schroder et al. 
(2003). The above equation is the best match to the mass-loss 
rates of Galactic disk C stars belonging to the Le Bertre (1997) 
sample with moderate mass loss (with colours in the range 2.5 5= 
J — K % 4), that is before the full development of the dust-driven 
superwind. 

In turn, the J - K colour for these models is computed with 
the aid of the relation 

(J - K) = 17.32 - 4.56 log T eB + 0.052 C/O (22) 

proposed by Marigo et al. (2003). The above equation extends 
the original colour- L e ff transformation for C stars by Bergeat et 
al. (2001), by introducing a C/O-dependence to account for the 
effect molecular blanketing on the spectra of these stars. 

Then, once the luminosity of a carbon-rich model overcomes 
Lsw the dust-driven superwind is assumed to start. This phase is 
described on the base of pulsating wind models presented by 
Wachter et al. (2002), in which the driving mechanism of mass 
loss is the radiation pressure on dust grains. From those models 
we derive a fitting relation, expressing the mass-loss rate as a 
function of stellar luminosity L, effective temperature T e s , mass 
M, pulsation period P, and photospheric C/O ratio (the piston 
amplitude is assumed Av p = 5 km s for all models) 

\o E [MI(M Q yx- 1 )] = -4.529 (23) 
-6.849 log(L eff /2600K)+ 1.527 log(L/10 4 L o ) 
-1.997 log(M/M ) - 0.995 log(P/650 days) 

+0 .6721og(£f ) . 

The numerical coefficients are quite similar to those of the 
fitting formula for M already proposed by Wachter et al. (2002; 
see their equation 1) obtained from the same set of wind models. 
We note that Eq. d23l contains an explicit dependence of M on 
the C/O ratio that is neglected in the Wachter et al. formula. Our 
choice is motivated by the fact that, though the influence of C/O 
is expected to be rather small for relatively low ratios (< 2 as is 
the case of Galactic C stars), it may instead become important 
for larger C/O ratios, e.g., in carbon stars belonging to metal- 
poor populations. 

A few examples of the predicted mass-loss rates are dis- 
played in Figs. Q3] and [16] In all cases, regardless of the C/O 
value, we note that i) the mass-loss rate is modulated by the 
variations of luminosity and effective temperature driven by the 
He-shell flashes, ii) most of the stellar mass is lost during the 
quiescent stages of high luminosity preceeding the occurrence 
of the thermal pulses, and iii) the super-wind regime is attained 
during the FM pulsation (P = Po). 

3. Basic model predictions 

The synthetic model just described in Sect. [2] has been ap- 
plied to calculate the TP-AGB evolution of an extended 
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set of stellar models, with initial masses in the inter- 
val 0.5 - 5.0 Mq, and for 7 choices of metallicity, Z = 
0.0001, 0.0004, 0.001, 0.004, 0.008, 0.019, 0.03. The TP- 
AGB phase is followed from the first thermal pulse up to the 
complete ejection of the envelope by stellar winds. 

In the following we describe the basic properties of our TP- 
AGB models, emphasizing the new aspects that derive from a 
number of updates applied to our synthetic code (described in 
Sect. |2j, the most important of which are: (i) an improved de- 
scription of the pulsation periods, so as to account for the switch- 
ing from the FOM to the FM of pulsation; (ii) a better recipe to 
predict the mass-loss rates, distinguishing between O-rich and 
C-rich models; (iii) the use of variable molecular opacities; (iv) 
a more detailed treatment of the third dredge-up and HBB. 

3.1. Transition luminosities 

The classical picture of TP-AGB evolution is one of luminosity 
generally increasing with time as the core mass grows along the 
evolution. This picture has already been demonstrated to fail un- 
der a series of circumstances (e.g., flash-driven luminosity vari- 
ations, very deep dredge-up, extinction of HBB, etc.), however 
it is still useful to define a few transition luminosities that mark 
the onset of particular regimes during the TP-AGB phase. The 
several panels in Fig. [T9l show a few characteristic luminosities 
as a function of the stellar mass, for the 7 metallicities explored 
in this work. These critical luminosities correspond to the quies- 
cent pre-flash values at the stages of: 

- The first termal pulse, L\jp\ 

- The switch from fundamental mode to first-overtone mode 
pulsation, Lp x - Lp ; 

- The transition from oxygen-rich to carbon-rich surface 
chemical composition, Lm-c; 

- The onset of the superwind phase of mass loss, Lsw, with 
M > 5 1O~ 7 M yr -1 ; 

- The luminosity at the end of the AGB phase, LAGB-tip- 

From comparing the various panels in Fig. [19] we derive the 
following indications: 

- The "chemical" transition from M- to C-type takes place at 
lower luminosities with decreasing Z. Moreover, the mini- 
mum mass to become a carbon star decreases with metal- 
licity, being ~ 2.0 M for Z = 0.019 and ~ 0.8 M for 
Z < 0.001. 

- The transition from the FOM to the FM pulsation occurs 
at higher luminosities with decreasing Z. In particular low- 
mass stars (say with M\ S; 1 - 1 .5 M ) of low Z enter the 
TP-AGB phase as FOM pulsators, while they should already 
appear as FM pulsators at solar or super-solar metallicity. 
At larger masses the transition to FM pulsation always takes 
place during the TP-AGB phase, regardless of Z. 

- Combining the above trends for Lp^p a (M\,Z) and 
Lm-c(M[,Z) it follows that at lower metallicities car- 
bon stars of any mass should first populate the FOM 
sequence, and then jump to the FM sequence in the P - L 
diagram. At increasing metallicities the Lm-c transition 
tends to occur at higher luminosities at given stellar mass, 
so that carbon stars should predominantly appear as FM 
pulsators. 

- The Lsw for both oxygen-rich and carbon-rich models is al- 
ways attained after the transition from FOM to FM pulsation, 
and it practically determines the final luminosity, Li P -agb at- 
tained at the termination of the AGB phase. 



- For carbon-rich models Lsw is very close to the Lm-c transi- 
tion at higher metallicities, e.g., as a star with initial Z > 

0. 019 becomes a carbon star the significant decrease in 

favours the quick development of the dusty superwind 
regime. At decreasing metallicities the surface cooling effect 
caused by the transition from M- to C-type domain is less 
pronounced, so that the onset of the superwind is delayed. 

- For oxygen-rich models Lsw takes place, on average, at 
larger luminosities with decreasing Z due to the higher T e g. 
We note that at very low metallicities, say Z < 0.001, low- 
mass models with M; < 0.8 M do not reach the minimum 
mass-loss rates (« 5 1O~ 7 M yr _1 ) required for the devel- 
opment of the dusty superwind. In these cases the reduction 
of the envelope mass and the consequent termination of the 
AGB phase are mostly determined by the core mass growth, 

1. e. the rate of outward displacement of the H-burning shell. 

3.2. TP-AGB tracks in the H-R diagram 

Figure [20] displays the complex morphology of the TP-AGB 
tracks in the H-R diagram at varying initial stellar mass and 
metallicity. Again, for the sake of simplicity, only quiescent 
stages - i.e. pre-flash luminosity maxima - are considered. 

As long as the mass loss remains so modest that the TP- 
AGB evolution can be safely considered at constant mass, the 
main parameter that affects the evolution in effective temperature 
is the surface chemical composition, specifically the C/O ratio 
(see Sect |3.3l >. Lower-mass stars that do not experience the third 
dredge-up, hence always remaining O-rich, are characterised by 
a slight decrease in L e ff as they evolve at increasing luminosi- 
ties along their tracks, with a slope lilogL/iilog L e (f essentially 
determined by the stellar mass and initial metallicity. 

This behaviour changes drastically as a model makes the 
transition from the O-rich to the C-rich configuration, which is is 
accompanied by a marked displacement toward lower effective 
temperatures (thick lines). As already discussed, this is the di- 
rect result of the enhanced molecular opacities as soon as as the 
photopsheric ratio passes from C/O < 1 to C/O > 1 . We also note 
that at decreasing metallicity, on average, the C-rich tracks de- 
scribe a wider excursion in T e g, compared to C-rich models with 
higher Z. This is due to the fact that at lower Z the transition 
to the C-star domain takes place earlier, so that models leave the 
O-rich branch toward lower L e ff s already during the initial stages 
of the TP-AGB. 

One remark should be made at this point. There are some 
intermediate mass models that, after having become carbon- 
rich as a consequence of the third dredge-up, are re-converted 
to the oxygen-rich configuration due to the prevailing effect 
of HBB (see Sect. |3.3l l. This is the case, for instance, of the 
(4M ,Z = 0.008) and (3.5M ,Z = 0.004) models. As shown 
in Fig.[20]the consequence of such sequence of chemical transi- 
tions (C/0< 1 — >C/0> 1 — >C/0 < 1) produces visible loops in 
the H-R diagram. 

At increasing mass and metallicity, say M\ > 4M and 
Z > 0.004, TP-AGB tracks do not describe a significant cool- 
ing in 7^ as HBB prevents the C/O ratio to increase above unity 
(except in the very last stages). However, we should remark that 
at very low metallicities (Z < 0.001) the role of HBB in keeping 
the surface C/O below 1 may not be effective any longer. Rather, 
efficient nuclear burning at the base of the convective envelope 
may favour itself the formation of C-rich models via the activa- 
tion of the ON cycle, as oxygen begins to be efficiently burnt in 
favour of nitrogen (see Fi g . [2TT> . Due to this, metal-poor tracks of 
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Fig. 19. Transition luminosities as a function of initial stellar mass and metallicity, 
namely: Ltp.i at the first thermal pulse; Lp^p at the switch from FOM to FM pul- 
sation; Lm-c at the transition from O-rich to C-rich surface chemical composition; Lsw 
at the onset of the dust-driven superwind regime; L max at the tip of the AGB. For the 
sake of simplicity, the plot is constructed considering the maximum quiescent lumi- 
nosities reached just before the occurrence of thermal pulses, even if the transitions 
may actually occur at some other phase of the corresponding pulse cycles 
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larger masses describe in Fig.[20]pronounced excursions toward 
lower r e ff. 

Finally, it is worth noticing that the luminosity evolution of 
more massive models is significantly affected by the occurrence 
of HBB. This process makes these models to be over-luminous 
compared to what predicted by the core mass-luminosity rela- 
tion in quiescent conditions. As we see in fact in Fig. [20] de- 
spite their core mass does not essentially grow because of the 
deep dregde-up, the corresponding tracks are characterised first 
by a relatively steep rise in luminosity as a consequence of HBB. 
Then, after reaching a maximum the tracks exhibit a final decline 
in luminosity, finally converging to the values predicted by the 
core mass-luminosity relation. This takes place at the onset of 
the superwind regime of mass loss, when the drastically reduced 
envelope mass determines the extinction of HBB. 



3.3. Surface chemical composition 

As already mentioned in Sect. 13.21 one key parameter that con- 
trols the morphology of the TP-AGB tracks in the H-R diagram 
is the photospheric C/O ratio, mainly affecting the evolution in 
T e g . This latter may remain unaltered during the whole TP-AGB 
phase - as is the case of lower mass models -, or it may be af- 



fected by the third dredge-up and HBB - as in the most massive 
models (refer to Sect. 12. 6b . 

It is interesting to compare a few emblematic cases, illus- 
trating the evolution of the CNO abundances and the C/O ratio 
in the convective envelope of models with different masses and 
metallicities. The results are schematically summarised as fol- 
lows. 

Models that experience only the third dredge-up show a pro- 
gressive increase of the I2 C abundance, hence of the C/O ratio 
which may overcome unity thus leading to the formation of car- 
bon stars. This case applies, for instance, to the M; = 2.5 M Q 
models shown in the bottom panels of Fig. [21] Note that the du- 
ration of the C-star phase is much longer for Z = 0.008 than for 
Z = 0.019. No important variations are expected in the abun- 
dances of N and O. 

In models that experience both the third dredge-up and HBB 
the evolution of the surface CNO abundances depends on the 
competition between the two processes. We distinguish three 
possible channels. 

1) If the third dredge-up initially dominates, then the C/O 
ratio is allowed to increase and the model enters the C-star do- 
main. The situation may be later reversed as HBB becomes more 
and more efficient in converting C in favour of N, so that the C/O 
ratio falls below unity. This case applies to the models with in- 
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Fig. 20. TP-AGB tracks in the H-R diagram as a function of stellar mass and metal- 
licity, as labelled. Luminosity and effective temperature correspond to the quiescent 
values just before the predicted occurrence of thermal pulses. Surface O-rich and C- 
rich configurations are distinguished with thin and thick lines, respectively. 



termediate mass, like the M\ = 3.5M , Z = 0.004 one shown in 
the middle-left panel of Fig. [21] 

It is worth remarking here that the evolution of these mod- 
els, hence their nucleosynthesis, strongly depends on the adopted 
molecular opacities and mass-loss rates. In fact, as shown by 
Marigo (2007), the use of a mass-loss formalism highly sensi- 
tive to T e g, e.g. Vassiliadis & Wood (1993), coupled to variable 
molecular opacities could weaken, or even quench, HBB shortly 
after the surface C/O ratio has overcome unity as a consequence 
of the third dredge-up. In that case, C/O does not decrease any 
longer and the model is expected to remain C-rich until the end 
of its evolution (see figure 2 in Marigo 2007). This specific point 
deserves more detailed investigation. 

2) The early transition to the C-star phase may be even pre- 
vented in models of larger mass when HBB is already strong 
enough to burn the newly dredged-up 12 C. These models are then 
predicted to keep the surface C/O ratio below unity for the whole 
TP-AGB phase, possibly except for the very last stages of in- 
tense mass loss when HBB stops and the third dredge-up makes 
them appear as obscured carbon stars (van Loon et al. 1997; 
Frost et al. 1998). This picture is exemplified, for instance, by 
the Mi = 4.5M , Z = 0.004 and Mi = 5.OM , Z = 0.008 mod- 
els displayed respectively in middle-right panel and bottom-left 
panel in Fig.l2Tl 



3) At very low metallicities HBB may even concur - together 
with the third dredge-up - to favour the formation of massive 
carbon stars. In fact, the high temperatures (^ 7.9 10 7 K) reached 
at the base of the convective envelope lead to the activation of the 
ON-cycle, with the consequent conversion of oxygen into nitro- 
gen. In this case the rise of the surface C/O ratio is mainly due 
to the depletion of the initial abundance of oxygen. According 
to our calculations this condition is met by models with M; <: 
4.0 M Q and Z < 0.001 (see Ventura et al. 2002 for similar re- 
sults). An example is given by the M\ = 5.0M Q , Z = 0.0004, 
shown in the bottom-right panel of Fig. |2T] We see that, after an 
initial decrease of carbon, the CN cycle attains its equilibrium 
regime - i.e. both 12 C/ 13 C and C/N reach their equilibrium val- 
ues - and on average carbon starts to increase steadily. Then, the 
activation of the ON cycle leads to a significant drop of oxygen 
which makes the C/O ratio increase above unity. Subsequently 
the ON-cycle goes to equilibrium as well so that carbon keeps 
more abundant than oxygen despite the steady increase of this 
latter. 



3.4. Duration of the TP-AGB phase 

Figure [22] shows the predicted TP-AGB lifetimes for both the 
oxygen-rich (with Mboi < -3.6) and C-rich phases, tm and re, 
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Fig. 21. Evolution of surface CNO abundances (in mole gr ') and C/O ratio during the TP-AGB phase of a few models with various 
initial masses and metallicities. 




Fig. 22. Predicted duration of the phases spent by TP-AGB models as M-type stars (tm: left panels) and C-type stars (jc- right 
panels), as a function of both stellar mass and metallicity. For the sake of comparison with observations we consider stellar models 
with luminosities higher than log L/L Q = 3.33, i.e., brighter than the RGB-tip at Mboi = -3.6 



as a function of the stellar mass, for all metallicities here consid- 
ered. It should be remarked that these results are derived from 
the calibration of our TP-AGB tracks, which will be discussed 
in Sect. |4] Moreover, the TP-AGB lifetimes displayed in Fig. [22] 
account also for the last evolutionary stages characterised by the 
superwind mass-loss rates, i.e. the contribution from the dust- 



enshrouded, optically-obscured evolution is included in the com- 
putation of T. 

The termination of the TP-AGB phase, hence its duration, is 
determined by all those factors that concur to make the star lose 
its envelope. Clearly the main agent is mass loss by stellar winds, 
(and its complex dependence on basic stellar parameters like 
luminosity, effective temperature, chemical composition, pulsa- 
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tion, dust formation, etc.), but a certain role is also played by the 
rate of outward advancement by the H-burning shell that reduces 
the envelope mass from inside (and its dependence on factors 
like the metal content, and the occurrence of the third dredge-up 
and HBB). Actually it is worth mentioning that for Z < 0.004 
low-mass TP-AGB models, say with M < 0.7M Q , never reach 
superwind mass-loss rates so that they leave the AGB because 
their envelope is eaten by the growth of the H-exhausted core. 

The behaviour with metallicity of tm at given mass, dis- 
played in the left panel of Fig. [22] depends on wether a TP-AGB 
model undergoes or not the third dredge-up. Lower-mass models 
that keep C/O < 1 up to the end of the AGB tend to have longer 
tm at decreasing metallicity because of the lower rates of mass 
loss and/or core-mass growth. 

In models of higher mass, that experience the third dredge- 
up, tm is crucially affected by this latter process - which deter- 
mines the stage at which the surface C/O > 1 - while it is es- 
sentially independent of the adopted mass-loss prescription. For 
these models the resulting trend of tm is controlled by the mass 
and metallicity dependence of the main two dredge-up parame- 
ters, A and M™ n . 

As we see there is a general tendency to have shorter tm 
at decreasing metallicity because of the earlier onset and larger 
efficiency of the dredge-up at lower Z. It is interesting to note 
that for Z <; 0.008 the tm(M0 curve displays a peak at about 
Mi <a 1 .9 - 2.2 Mq. This finding is to be at least partly ascribed to 
the relation M c .itp(M) at the first thermal pulse, which presents a 
minimum just in correspondence to these stellar masses, so that 
a relatively large growth of the core mass is required before the 
M™ n is attained. As we discuss in Sect. 14.21 this feature is in 
agreement with the observed maximum of M-type star counts in 
the Magellanic Clouds' clusters. 

Concerning tc, shown in the right panel of Fig. [22] we note 
that the mean trend with metallicity is reversed compared to 
tm, i.e. the C-rich phase is expected to last longer at decreasing 
metallicity. For instance at M\ = 2.5 M , tc passes from 0.6 Myr 
for Z = 0.019 to 2.8 Myr for Z = 0.008. However, we also 
see that moving to very low Z, tq does not increase drastically, 
rather it keeps on average almost invariant, e.g., tc ~ 2.4 Myr 
for Z = 0.0004. This is the net result of the interplay between 
mass loss and dredge-up through their dependence on metallic- 
ity. 

3.5. Mass loss rates, periods and luminosities of variable 
AGB stars 

The use of theoretically-derived mass loss rates constitutes a ma- 
jor change with respect to our previous works. Leaving a detailed 
discussion of this subject to a future paper, below we simply il- 
lustrate the implications of our present choice. 

Figure [23] compares a few of our O-rich TP-AGB tracks of 
solar composition with the data for Galactic M-type Miras in the 
M-P diagram. It can be noticed that the models reproduce well 
the ranges of mass-loss rates and periods covered by the data. In 
particular, the significant dispersion in M at given P exhibited 
by the Miras is accounted by the models at varying mass. The 
observed spread cannot be obtained, instead, by applying, for 
instance, the Vassiliadis & Wood (1993) formalism that gives 
M as a sole function of P. Moreover, the onset of the superwind 
phase (with M a 10~ 5 -10 -4 M yr _1 ) is naturally predicted with 
the adopted scheme and allows the reproduction of the OH/IR 
data over a wide period range. 

Similar comments can be made about mass loss rates on 
the C-rich phase, as illustrated in Fig. [24] What is more evi- 
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Fig. 23. Mass-loss rate versus period for O-rich Mira variables. 
Observed data (circles) correspond to Galactic M-type Miras (Le 
Bertre & Winters 1998; Groenewegen et al. 1999). For the sake 
of comparison we show the semi-empirical relation for Miras 
proposed by Vassiliadis & Wood (1993; dashed line), together 
with the theoretical predictions for evolving O-rich TP-AGB 
models with Z = 0.019 and various initial mass, based on Bowen 
& Willson (1991) wind models. 
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Fig. 24. Mass-loss rate versus period for C-rich Mira vari- 
ables. Observed data correspond to Galactic C-rich Miras from 
Groenewegen et al. (1998; squares), and Schoier & Olofsson 
(2001; triangles). For the sake of comparison we show the semi- 
empirical relation for Miras proposed by Vassiliadis & Wood 
(1993; dashed line), together with the theoretical predictions for 
evolving C-rich TP-AGB models with Z = 0.019 and various 
initial mass, based on Wachter et al. (2002) wind models. 



dent in this case is that models reproduce well the plateau at 
M a 10" 5 M o yr _1 in M, corresponding to the superwind phase 
of C-type stars. 

Extensive observations of Miras in the LMC have shown that 
these variables populate a well-defined period-luminosity (P-L) 
relation (e.g., Feast et al. 1989; Hughes & Wood 1990; MACHO 
survey), as displayed in Fig. [25]for both oxygen- and carbon-rich 
objects. 
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Fig. 25. Period-luminosity rela- 
tion of Mira variables in the 
LMC. Observed data correspond 
to M-type (empty circles) and C- 
type Miras (filled circles) LMC 
(Feast et al. 1989; Hughes & 
Wood 1990), assuming a dis- 
tance modulus fio = 18.5. Lines 
show the theoretical tracks de- 
scribed by TP-AGB models for 
two choices of metallicity and 
with initial stellar masses rang- 
ing from 0.6 to 5.0 M G (a few 
values are indicated). It is as- 
sumed that pulsation occurs ini- 
tially in the FOM (P\ sequence) 
and then it switches to the FM 
(Pq sequence). See text for more 
esplanation. 



TP-AGB evolutionary tracks at varying initial mass and with 
Z = 0.019 and Z = 0.008 are superimposed for the sake of 
comparison. According to the assumed description of pulsation 
(see Sect. 12. 7K tracks split into two bundles in the P — L diagram, 
the one at shorter P corresponding to the FOM, and the other 
at longer P corresponding to the FM. Clearly this latter group 
turns out to be consistent with the observed location of Miras, 
supporting the indication that these variables are pulsating in the 
fundamental mode. 

A detailed investigation of the several variability sequences, 
in addition to that of Miras, that populate the P—L diagram (e.g. 
Wood et al. 1999) are beyond the scope of the present study. To 
this aim it is necessary to perform population synthesis simula- 
tions that couple updated TP-AGB tracks with LPV models for 
different pulsation modes. This kind of analysis is underway. 

3.6. The initial-final mass relation 

The initial-final mass relation (IFMR) is affected by the compe- 
tition between core growth by shell burning and core reduction 
by third dredge-up events, the final result depending on the ter- 
mination of the AGB evolution by stellar winds. Since dredge- 
up is more efficient at higher masses and lower metallicities, one 
may naively expect that this competition will favour flatter IFMR 
at lower metallicities. Actually, this may not happen because the 
IFMR also carries a strong imprinting of the core mass at the 
first thermal pulse, which is in general larger at lower metallic- 
ities (see e.g., figure 1 in Girardi & Bertelli 1998). Moreover, 
the larger efficiency of HBB at lower metallicities may reduce 
the TP-AGB lifetime in such a way that less time is left for the 
core mass to significantly change along the TP-AGB. Finally, 
the IFMR is modulated by the TP-AGB lifetimes, the complex 
dependence of which as a function of mass and metallicity is 
shown in Fig. [22] Therefore, the actual shape of the IFMR and 
its variation with metallicity depends on several parameters and 
it cannot be easily predicted via simple reasonings. 

Figure [26] shows the IFMR as predicted by present models 
of various metallicities, compared to recent observational data 
for white dwarfs in binary systems and open clusters. In gen- 
eral, both the cluster ages and white dwarf progenitor masses 
have been estimated (or checked) by using tracks and isochrones 
from Girardi et al. (2000), so that the plot should not contain in- 
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Fig. 26. The initial-final mass relation from our models (lines) 
for several metallicities, together with the empirical data for 
white dwarfs in open clusters and binary systems: the Hyades 
(full circles, Claver et al. 2001), M 37 (open squares, Kalirai et 
al. 2005), Praesepe (crosses, Claver et al. 2001, and Dobbie et al. 
2006), M 35 (full squares, Williams et al. 2004), Sirius B (star, 
Liebert et al. 2005a), Pleiades (full triangle, Dobbie et al. 2006), 
and NGC 2516 (open circles, Ferrario et al. 2005). The gray area 
indicates the position (mean and lcr) of the primary mass peak 
of field white dwarfs as measured by Liebert et al. (2005b), that 
we indicatively associated to the 1.0-1.5 M interval. 



consistencies associated to the use of different prescriptions for 
convective overshooting. It is worth noticing that: 

1) The predicted IFMR have a significant dependence on 
metallicity. For a given initial mass the allowed range of final 
masses in about 0.1 M e wide. The dependence is in the sense 
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the IFMR is steeper at low metallicities. Although the final mass 
generally increases with the initial mass, there are limited mass 
ranges in which this trend is reversed (e.g., for the Z = 0.03 
tracks in the 2.2-3 M interval), as well as small track-to-track 
variations that arise from the complex interplay between the evo- 
lutionary parameters involved. 

2) The data seems to be too dispersed to allow a firm testing 
of the model predictions, but this dispersion largely depends on 
the data for the M 37 cluster (from Kalirai et al. 2005), which - 
apart from a few low-mass outlyiers - lie systematically above 
the mean relation drawn by the remaining objects. As remarked 
by Dobbie et al. (2006), the mass determinations of M 37 white 
dwarfs may still change considerably since its initial metallicity 
is not well determined. Excluding the M 37 white dwarfs, the 
IMFR is well defined for the complete range of M\ > 2.5 M Q . 
On the other hand, good-quality data covering the 1.5 - 2.5 M Q 
mass range are presently not available. Notice that we have rep- 
resented the < 1 .5 M Q mass range by field white dwarfs, whose 
initial mass is not directly measured. 

3) If we disconsider the M 37 white dwarfs, the comparison 
between the data and the model IFMRs for close-to-solar metal- 
licities (the Z = 0.019 and 0.03 curves) turns out to be very satis- 
factory over the entire M; < 5 M G mass range. It seems however 
advisable to extend the TP-AGB model computations to slightly 
larger initial masses, in order to explain the few points observed 
at Mi > 5.5 M Q . In this regard, more recent evolutionary tracks 
from Padova (Bertelli et al., in preparation) suggest that the high- 
est mass for entering in the TP-AGB phase, using Girardi et al.'s 
(2000) prescription for overshooting, is located somewhere be- 
tween 5.5 and 6 M Q , rather than equal to the 5 M Q mentioned by 
Girardi et al. 

Future works in this series will improve on this latter point, 
and present detailed comparisons with the mass distribution of 
field white dwarfs. 



4. Calibration of TP-AGB model parameters 

In this section, we first briefly present the set of Magellanic 
Cloud data used to constrain the TP-AGB model parameters, to- 
gether with the methods employed to simulate the data starting 
from the TP-AGB tracks. The calibration process is fully de- 
scribed in the Appendix. 

We remark that previous sets of calibrated TP-AGB models 
(e.g., Groenewegen & de Jong 1993; Marigo et al. 1999) have 
used the C-star luminosity functions in the Magellanic Clusters 
as, basically, the only observational constraint. To the CSLFs, 
we now add the C- and M-type lifetimes as to be described be- 
low. Previous works have also adopted too simple descriptions 
for the star formation rate (SFR) and age-metallicity relations 
(AMR) of both galaxies, like for instance a constant SFR, and 
a constant metallicity with varying age. We now improve upon 
this approach by adopting SFR and AMR independently derived 
by other authors to simulate the data. 

4.1. Carbon star luminosity functions in the Magellanic 
Clouds 

The CSLF in several Local Group galaxies have been re- 
viewed and intercompared by Groenewegen (2002; see also 
Groenewegen 1999), who kindly provided us with his latest data 
files for both Magellanic Clouds. The observed CSLF in the 
LMC makes use of the Costa & Frogel (1996) m\, \ data, which 
refers to 895 carbon stars among the original sample of 1035 



stars distributed in 52 LMC fields over an area of 8° x 8°. This 
sample of carbon stars can be considered almost complete and 
not limited by apparent magnitude, since the faintest carbon stars 
are at least 2 mag brighter (I ~ 15) than the faint limit of detec- 
tion (7 ~ 17) of the original surveys (Blanco et al. 1980, Blanco 
& McCarthy 1983). The fraction of optically-obscured carbon 
stars, which would have been missed in the survey because of 
dust obscuration, has been estimated to be lower than 3 % for 
Mroi < -6 (Groenewegen &de Jong 1993). The conversion from 
the photometry to bolometric magnitudes has been performed by 
Costa & Frogel (1996), who mention a mean error on the eval- 
uation of mboi of approximately 0.34 mag. Finally, the absolute 
bolometric magnitudes, Mroi, are then obtained by adopting a 
distance modulus of fj.o = 18.5 for the LMC. The resulting CSLF 
extends approximately from M\, Q \ = -3 up to M\, D i = -6.5, with 
the peak located at around M\, D i = -4.875. 

The observed CSLF in the SMC has been derived by 
Groenewegen (1997) from a sample of 1636 stars observed by 
Rebeirot et al. (1993), and adopting bolometric correction by 
Westerlund et al. (1986). The adopted distance modulus for the 
SMC is //o = 19.0. The first bin of this CSLF contains all stars 
fainter than Mx>i = -2.5, which comprise less than 3 % of the 
total sample of SMC carbon stars. Their intrinsic magnitudes 
can be as low as Mx>i = -1.8, as found by Westerlund et al. 
(1995). The formation of carbon stars at such low luminosities 
can hardly be understood as the result of the third dredge-up in 
single AGB stars, and most likely they arise from the transfer of 
mass in close binary systems (see e.g. de Kool & Green 1995; 
Frantsman 1997). These low-luminosity carbon stars will not be 
considered in this work, since we are dealing with single star 
evolution. For all practical comparisons, we will consider the 
CSLF in the SMC as extending down to Mroi — -2.5. The re- 
sulting CSLF extends up to Mboi = -6.5, with the peak located 
at Mboi = -4.375. 

Recently Guandalini et al. (2006) have claimed that bolo- 
metric corrections for C stars are systematically underestimed 
because neglecting the flux in the mid-IR, and suggested that ac- 
counting for this flux would imply overall much brighter CSLFs. 
Although their arguments are unimpeachable, it is worth notic- 
ing that their study is based on a sample of C stars in the Milky 
Way selected to have mid-IR photometry. Such a sample is much 
more prone to contain IR-bright objects - of cool temperatures 
and undergoing significant mass loss and dust obscuration -, 
than samples extracted from optical to near-IR surveys. This is 
also indicated by the fact that the majority of their stars have 
J - K > 2.0 (see their figure 7), which corresponds to the colour 
interval of C stars strongly affected by circumstellar dust. The 
bulk of C stars observed in the Magellanic Clouds, instead, have 
J-K comprised between 1 .2 and 2.0 (e.g. Nikolaev & Weimberg 
2000), and it is expected to be overall much hotter and cleaner 
from dust than Guandalini et al.'s (2006) sample. Moreover, we 
remark that both observations (van Loon et al. 2005) and popu- 
lation synthesis models based on the present evolutionary tracks 
(Marigo et al., in prep.) agree in estimating as « 10 - 20% the 
fraction of dust-enshrouded C-stars in the Magellanic Clouds, 
i.e. the bulk of their C-star population consists of optically- 
visible objects with just mild mid-IR emission. Therefore, we 
consider it unlikely that the consideration of the mid-IR flux 
would significantly affect the CSLFs here discussed. 

4.2. TP-AGB lifetimes: M- and C-type phases 

Limits to the lifetimes of the M- and C-type TP-AGB phase have 
been recently derived by Girardi & Marigo (2007), using avail- 
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able data for C and M giants with Mboi < -3.6 in Magellanic 
Cloud clusters (Frogel et al. 1990), together with the clusters' 
integrated V-band luminosities (from Bica et al. 1996 and van 
den Bergh 1981). They essentially find the C-star phase to have 
a duration between 2 and 3 Myr for stars in the mass range from 
~ 1.5 to 2.8 M , with some indication that the peak of C-star 
lifetime shifts to lower masses as we move from LMC to SMC 
metallicities. The M-giant lifetimes is found to peak at ~ 2 M 
in the LMC, with a maximum value of about 4 Myr, whereas 
in the SMC their lifetimes are poorly constrained by the data. 
Therefore, in the following we will only deal with the C-star 
lifetimes for both LMC and SMC, and with the M-star lifetime 
above M bo i < -3.6 for the LMC. 

Our aim is to obtain theoretical TP-AGB lifetimes that fit 
the Girardi & Marigo (2007) values, within the 67 % confidence 
level of the observations (lcr for the most populated age bins), 
over the complete 0.8 - 5 M interval for the formation of TP- 
AGB stars. Although this 67 % confidence level interval may 
look quite ample, they constitute very useful constraints. In fact, 
Girardi & Marigo (2007) show that many TP-AGB models in 
the literature do present too short C-star lifetimes, and/or fail to 
reproduce the observed trends of the lifetimes with the initial 
mass. 

4.3. Modelling the CSLF and lifetimes 

In order to simulate the CSLF starting from TP-AGB tracks, we 
follow a complex but well-tested procedure, which is somewhat 
more elaborated than the one described in Marigo et al. (1999). 
First, the tracks on their quiescent phases of evolution, are as- 
sembled together with the previous evolution given by Girardi et 
al. (2000) tracks. Theoretical isochrones are then constructed via 
interpolation for any intermediate value of age and metallicity 
(Marigo & Girardi, in preparation). The isochrone section corre- 
sponding to the C-star phase - with their quiescent luminosities 
- is then isolated; for each point of it, the pulse-cycle luminosity 
variations are reconstructed using Wagenhuber & Groenewegen 
(1998) formula together with basic stellar parameters like the 
metallicity and core mass. This allows us to attribute to each 
quiescent point in the isochrone the complete luminosity proba- 
bility distribution function, and hence to construct the CSLF for 
each single-burst stellar population via integration along the cor- 
responding isochrone. Then, we integrate over population age, 
weighting each age by its relative SFR and using the AMR to 
select the right metallicity value at each age. 

As for the LMC, we have adopted the SFR determined by 
Holtzman et al. (1999) in an outer field close to Hodge 4, us- 
ing deep HST data. For the SMC, we have used the global SFR 
as determined by Harris & Zaritsky (2004) from the Magellanic 
Clouds Photometric Survey. In both cases, the SFR have been re- 
constructed using objective algorithms of CMD fitting. Their re- 
sults are mostly determined by the star counts in the main phases 
of pre-AGB evolution (main sequence, subgiant and red giant 
branches, red clump) and, importantly to us, are unlikely to be 
sensitive to the (relatively very few) TP-AGB stars present in 
their data. 

We adopt the AMRs from the "bursting models" of chem- 
ical evolution for the LMC and SMC by Pagel & Tautvaisiene 
(1998). These relations, although of theoretical nature, are de- 
signed to reproduce the age-metallicity data for star clusters in 
both Clouds. 

The TP-AGB lifetimes in the C- and M-type phases are de- 
rived directly from the tracks, first identifying the turn-off mass 
and population age of each TP-AGB track, then interpolating 



the turn-off mass and TP-AGB lifetimes for each metallicity- 
age value given in the AMR. As a result, we derive the expected 
TP-AGB lifetimes (above M\, \ = -3.6) as a function of turn- 
off mass, that can be directly compared to the data tabulated by 
Girardi & Marigo (2007). 

We notice that, with the present prescriptions for the mass 
loss rates, the superwind phase (with M S 5 .6 x 1 0~ 7 M yr~ 1 ) for 
which the stars are likely to be optically obscured, lasts less than 
20 percent of the TP-AGB lifetime for all tracks with Z < 0.008. 
For models with 1.5 ^ (M/M ) £ 2.5 of such metallicities, 
which enter the superwind phase as C stars and are the most 
relevant for the comparison with MC data, the superwind phase 
lasts typically less than 15 percent of the C-star lifetime. These 
numbers are in rough agreement with those derived from ob- 
servations in the MCs (van Loon et al. 2005). A detailed anal- 
ysis of this point is left to subsequent papers in this series. We 
advance however that the neglection of objects that are likely 
optically-obscured in the theoretical CSLFs and lifetimes, as we 
verified, would not significantly impact on our calibration iter. 
The Appendix summarizes the calibration iter, starting from the 
past calibration of Marigo et al. (1999) and arriving at the final 
result depicted in Fig. [27] 

4.4. Implications from the calibration 

The calibration sequence described in the Appendix yields a 
number of interesting and useful indications, namely: 

1) The combined use of both CSLFs and stellar lifetimes in 
the Magellanic Clouds as calibrating observables strengthens the 
reliability of the underlying TP-AGB description. In fact, we find 
that sets of models, consistent with the data for the C- and M-star 
lifetimes, may instead not reproduce the corresponding luminos- 
ity functions and viceversa. 

2) Due to the pronounced metallicity dependence of the third 
dredge-up, any theoretical investigation of the C-stars' popula- 
tions should account for the history of chemical enrichment in 
the host galaxies, then relaxing the too simple assumption of 
constant metallicity. Indeed, the simultaneous reproduction of 
the CSLFs in both Magellanic Clouds provides a strong indica- 
tion in this respect, as we see in Fig. 127] when comparing the 
results for Z = const, to those obtained including proper AMRs 
in the simulated galaxies. In particular, we see that the excess 
of bright C stars and the deficiency of fainter objects, predicted 
while assuming Z = const., are removed when accounting for 
the metallicity-dependence of the third dredge-up. This means 
that the extent of the long-standing theoretical difficulty in pro- 
ducing faint carbon stars, known as 'the carbon star mystery" 
since Iben (1981), could be actually less important than so far 
believed, just because of a metallicity mismatch in the compar- 
ison between observations and predictions. Moreover, the fact 
that the populations of C stars in the two galaxies share a com- 
mon metallicity range, say around Z = 0.004, should make our 
metallicity calibration for the third dredge-up more constrained. 

3) The peak of the M-star lifetimes in the LMC clusters (lim- 
iting to stars brighter than the RGB-tip) is recovered only assum- 
ing that the third dredge-up in stars with M\ ~ 1.7 — 2.5M and 
Z a 0.008 does not take place as soon as they enter the TP-AGB 
phase, but somewhat later. The M™ n (M, Z) formalism derived 
by K02 turns out to be consistent with this indication. In this 
way, we alleviate the problem of under-estimating the M-star 
populations that affects the simulations presented in Cioni et al. 
(2006ab). 

4) The decreasing lifetimes of C stars with Mi > 2.5M seen 
in both Magellanic Clouds is due to the concomitant effects of 
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= f(M,Z) A = f(M,Z) P = P 1 -P 





Fig. 27. The final model calibration. Left panels: The observed C star luminosity functions (CSLF) in the Magellanic Clouds 
(Groenewegen 2002; shaded areas), as compared to the synthetic ones derived from present models. Right panels: The observed C- 
and M-type lifetimes in the Magellanic Clouds (Girardi & Marigo 2007), as compared to the same models. 



their brighter luminosities - that favour the earlier onset of the 
superwind - and, in models of larger masses (with M\ ^ 4.0M o ), 
of the occurrence of HBB - that may later reconvert a C star into 
an M star, or even prevent the transition to the C-star domain. 



5. Final remarks 

The set of TP-AGB evolutionary tracks presented in this paper 
constitutes a very useful database for population studies, due to 
several distiguishing features: 

1) We have used updated ingredients and improved treat- 
ment of some key physical processes in the synthetic AGB code, 
as extensively discussed in Sects. [2] and [3] The most impor- 
tant updating is no doubt the use of variable molecular opaci- 
ties instead of the scaled-solar tables still used in complete AGB 
models. Previous sets of models using these opacities have al- 
ready been used (Marigo 2002; Marigo et al. 2003; Cioni et 
al. 2006ab) and distributed by us (see Cioni et al. 2006a and 
http : //pleiadi . oapd . inaf . it), but this is the first time cal- 
ibrated TP-AGB tracks using this approach are released. 

2) The present model calibration (Sect. |4]i is based not only 
on the reproduction of the CSLFs in both Magellanic Clouds, but 
also on the data for C- and M-star counts in Magellanic Cloud 
clusters. It means that present models have not only the right 
luminosities of the TP-AGB phase (as ensured by the CSLF cal- 
ibration), but also that they present, as far as possible, the right 
lifetimes at Magellanic Clouds' metallicities. The contribution 
of any post-MS evolutionary stage to the integrated light of a 
single-burst stellar population is proportional to the product of 
luminosity and lifetime - also known as the nuclear fuel in the 
context of the fuel consumption theorem (Renzini & Buzzoni 
1986). Therefore, evolutionary population synthesis models that 
use the present TP-AGB tracks will present the right contribu- 
tion of TP-AGB stars to the integrated light at Magellanic Clouds 
metallicities, and hopefully, also a good description of the way 
this contribution varies with metallicity. 



3) The present tracks nicely complement the sets of low-and 
intermediate-mass models of Girardi et al. (2000) and Girardi 
(2002, unpublished), that deal with the evolution from the zero- 
age main sequence up to the first thermal pulse. In fact, when 
coupling both sets of tracks together one will find no discon- 
tinuities in quantities like the core mass and chemical surface 
composition. Discontinuities in the HR diagram may appear due 
to the use of slightly different low-temperature opacities; we 
checked however that they are tiny and certainly much less seri- 
ous than those one have when amalgamating tracks from differ- 
ent sources. Therefore, one has the possibility of building consis- 
tent tracks all the way from the ZAMS to the end of the TP-AGB. 
Of course, this sequence can also be extended with any suitable 
set of post-AGB tracks (e.g., Vassiliadis & Wood 1994), so as to 
cover the subsequent evolution up to the WD stages. 

4) We make a first attemp to describe the switching of pul- 
sation modes between first overtone and fundamental one along 
the TP-AGB evolution. These changes have some consequence 
in terms of the history of mass loss, and open new possibilities 
to better calibrate TP-AGB tracks by means of comparison with 
the extensive variability surveys in the Magellanic Clouds. 

Forthcoming papers will present extensive and ready-to-use 
theoretical isochrones and chemical yields from these tracks, as 
well as several population synthesis applications. 
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Appendix A: The model calibration 

This section summarises the calibration iter, mainly involving 
the parameterised decription of the third dredge-up. We also dis- 
cuss the strong influence played by other parameters, like the 
assumed pulsation mode on the mass-loss rates, and the age- 
metallicity relation of the host galaxy. 

The most relevant calibrating steps are illustrated through the 
sequence of panels from Figs. IA. ll to l27l where labels indicate 
the main model assumptions, namely: 

- Molecular opacities, either fixed for solar-scaled mixtures, 
Ktix, or variable chemical composition, /c var ; 

- Dredge-up efficiency A; 

- Onset of the third dredge-up, determined by either the tem- 
perature 7^ red parameter, or the classical M™ n parameter. 

- Pulsation mode, either fundamental with P = Po, or first 
overtone with P - Pi ; 

- Metallicity of the simulated galaxy, either assumed constant 
over its entire history, Z = const, or made vary according to 
a specified age-metallicity relation, AMR. 

Figure I A. 1 1 illustrates the results obtained with the same 
prescriptions for the third dredge-up as in our previous study 
(Marigo et al. 1999). The most significant difference between 
the two sets of models is the adoption in the present calcula- 
tions of variable molecular opacities in place of those for solar- 
scaled mixtures. In both cases mass loss is treated with the aid 
the Vassiliadis & Wood (1993) formalism, always assuming FM 
pulsation (P = Po). 

We notice that the combination (log T dled = 6.40, A = 0.5) 
now fails to reproduce the CSLF in the LMC - instead provid- 
ing the best fit in Marigo et al. (1999) previous calibration - 
while it well recovers the C-star lifetimes in both Magellanic 
Clouds. A reasonable fit to the CSLF in the LMC is obtained 
with the (log T dKd = 6.45, A = 0.5) case, but the duration of the 
C-star phase as a function of M\ is somewhat underestimated. 
The (log T dred = 6.50, A = 0.5) set of models yields poor re- 
sults in respect to all observables. Finally, we should remark that 
all combinations of dredge-up parameters do not reproduce the 
CSLF in the SMC, and they substantially under-estimate the du- 
ration of the M-star phase in the 1.8 - 3.OM range. 

Figure IA.2I is based on models calculated with exactly the 
same input parameters of Fig. IA.11 except that pulsation is as- 
sumed to occur in the FOM. With the assumption P = P\, the 
predicted CSLFs tend to populate brighter luminosities, as a con- 
sequence of the fact that shorter periods correspond to lower 
mass-loss rates, hence longer TP-AGB lifetimes. From the com- 
parison with the observed data it turns out that (i) none of the 
(T dled , A) combinations leads to a satisfactory reproduction of 
the CSLFs; (ii) the C-star lifetimes are relatively well recovered; 
(iii) the M-star lifetimes in the LMC are still notably underesti- 
mated in the 1.8 - 3.OM range. 

Figure IA.3I introduces a few important changes, namely (i) 
the assumption of constant A is relaxed to account for a more 
realistic dependence on mass and metallicity according to K02 
recipe; (ii) pulsation periods are calculated following the scheme 
indicated in Sect. 12.71 that predicts the mode switching from 
FOM (P = Pi) to FM (P = P ); and (iii) the inclusion of an 
age-metallicity relation in the galaxy simulations instead of con- 
stant metallicity. As a matter of fact, the use of a A{M, Z) law sig- 
nificantly improves the results for the CSLF in the LMC, while 
a sizeable overproduction of luminous C stars affects the CSLF 
in the SMC. Moreover, the problem with the too short M-star 
lifetimes still remains. 



The way out to this latter discrepancy is obtained when we 
adopt the Mf m (M,Z) formalism suggested by K02, in place of 
the temperature parameter T dred to determine the onset of the 
third dredge-up. This effect is illustrated in Fig. |A.4| The reason 
is that while the T dled criterion predicts that models with M\ <: 
2M Q experience the mixing events already from the first thermal 
pulse, with the M™ n (M, Z) scheme the onset of the dredge-up is 
delayed, then allowing for a longer duration of the M-type phase 
for these models. 

We should also note that the theoretical CSLFs obtained 
with Z = const, notably differ from the observed ones in both 
Magellanic Clouds. The results gets better when assuming a 
proper AMR for the host galaxy, but the comparison with the 
observed data is still not completely satisfactory. 

Finally, the best fit to the CSLFs, the M-type, and the C-type 
lifetimes in the Magellanic Clouds is reached with the aid of 
the scheme for the third dredge-up presented in Sect. 12.51 The 
original K02 formalism is partly modified in order to (i) correct 
for convective overshooting, (ii) anticipate the occurrence of the 
third dredge-up, (iii) and make it more efficient at lower metal- 
licities. The results are displayed in Fig. [27] 
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K var ^ = 0.5 P = P Z = const. 




Fig.A.l. Left panels: The observed C star luminosity functions (CSLF) in the Magellanic Clouds (Groenewegen 2002; shaded 
areas), as compared to the synthetic ones derived from models with different dredge-up parameters, assuming Vassiliadis & Wood's 
(1993) mass loss prescription and pulsation on the FM. Right panels: The observed C- and M-type lifetimes in the Magellanic 
Clouds (Girardi & Marigo 2007), as compared to the same set of models. 



K var A = 0.5 P = P, Z = const. 
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Fig. A.2. The same as Fig. lA. 11 but for models assuming pulsation on the FOM. 
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Fig. A.4. The same as Fig. IA.3l but for models assuming the M™ n (M, Z) formalism suggested by K02. 



